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\ Echocardiography  is  used  routinely  in  cardiology. 
The  typical  dynamic  range  encountered  spans  about  90  dB. 
The  objective  of  this  research  is  to  devise  a  means  to 
compress  this  wide  dynamic  range  to  a  narrower  dynamic 
range  for  grey  scale  display  or  for  digital  data  storage. 
An  adaptive  dynamic  range  algorithm  to  maximize  the 
information  content  of  the  echo  signal  is  developed  and 
tested  experimentally  and  also  using  simulations. 

An  echo  signal  model  is  developed  to  formulate 
the  dynamic  range  compression  algorithm  and  to  facilitate 
the  evaluation  of  dynamic  range  compression  algorithms 
with  simulated  data.  The  echocardiographic  signal  is 
modeled  as  a  nonstationary  stochastic  process  consisting 
of  specular  and  scatterer  components.  The  separation  of 
scatterer  components  from  blood  volume  is  important  in 
establishing  a  significance  signal  level.  A  lower  bound 
on  this  significance  level  based  on  acoustic  wave 
divergence  is  suggested.  The  significant  echo  signal 
of  specular  and  scatterer  components  possess  wide  dynamic 
range  throughout  the  scanning  volume  due  to  their  non- 
stationarity .  Based  on  this  model,  a  general  adaptive 

dynamic  range  compression  algorithm  is  developed.  Vs 
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The  adaptive  dynamic  range  compression  algorithm 
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relies  on  the  fact  that  the  stochastic  process  is  of  \the 
slowly  fluctuating  type.  Based  on  an  independent  estimate 


of  local  power  statistics,  the  algorithm  tends  to  normalize 
the  echo  signals  so  that  they  approximate  pseudo-stationary 
processes.  An  integral  part  of  the  algorithm  is  a  non¬ 
linearity  introduced  so  that  only  signals  of  significance 
are  normalized.  This  general  adaptive  dynamic  range 
compression  algorithm  also  applies  to  multidimensional 
echocardiographic  imaging. 

Two  specific  implementations  of  this  general  adaptive 
algorithm  are  suggested  for  M-mode  and  three-dimensional 
sector-scanned  snapshot  mode  echocardiographic  imaging 
systems.  Their  ability  to  adapt  to  the  nonstationary 
signals  and  to  suppress  insignificant  signals  is  examined. 
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CHAPTER  I 


INTRODUCTION 

Echocardiography,  a  pulse-echo  ultrasonic  technique, 
is  fast  gaining  favor  as  a  diagnostic  medical  instrument 
in  the  evaluation  and  analysis  of  heart  disease.  It  is 
comparatively  low  risk  due  to  its  nonionizing  and  non- 
invasive  nature.  Many  echocardiographic  systems  and 
techniques  have  been  developed  since  the  initial  use  of 
ultrasound  to  ascertain  the  volume  of  the  human  heart  in 
the  1950's.  This  thesis  considers  a  specific  problem 
encountered  in  echocardiography  —  that  of  dynamic  range 
compression.  A  general  dynamic  range  compression  technique 
applicable  to  multidimensional  (which  includes  time) 
echocardiographic  imaging  systems  is  presented.  The 
specific  application  of  such  a  technique  to  a  three- 
dimensional  snapshot  camera  is  emulated.  This  technique 
is  evaluated  with  simulated  echo  data  for  a  model  presented 
here.  Experimental  data  from  volunteers  are  also  presented 
for  comparison  with  the  model  and  evaluation  of  the  dynamic 
range  compression  algorithm. 

1.1  Historical  Remarks 

The  utilization  of  pulse  echo  ultrasonic  systems 
for  cardiology  started  in  the  early  1950's.  As  in  all 
ultrasonic  instruments,  the  evolution  was  stimulated  by  the 
development  of  sonar  during  World  War  II.  Firestone , 


in  the  early  1950 's,  started  the  use  of  pulse-echo  ultra- 

(24) 

sound  for  nondestructive  testing.  Keidelv  '  was  one  of 

the  first  investigators  to  use  ultrasound  to  examine  the 

heart.  His  technique  was  to  transmit  ultrasonic  waves 

through  the  heart  and  record  the  effects  of  the  ultrasound 

on  the  otherside  of  the  chest.  The  purpose  of  his  work  was 

to  try  to  determine  cardiac  volume. 

The  true  beginning  of  echocardiography  did  not 

(9) 

start  till  Hertz  and  Elder'  7  modified  a  commerical  non¬ 
destructive  testing  ultrasonoscope  in  1953.  Preceding 
and  during  this  period,  ultrasonic  diagnosis  also  flourished 
in  other  medical  applications.  One  of  the  dominant  figures 
in  diagnostic  ultrasound  was  D.  H.  Howry. A  period  of 
investigation  and  interpretation  of  the  echograms  followed 
Hertz  and  Elder’s  first  venture  into  echocardiography. 

The  display  utilized  at  that  time  was  the  A-mode  display. 

The  A-mode  display  presents  echo  strength  versus  range 
(actually  time)  data.  Since  the  pulses  were  repeated,  an 
additional  dimension  is  time.  Thus,  the  early  echocardio- 
graphic  imaging  systems  presented  the  echo  information 
along  two  independent  variables,  range  and  time. 

"Technical  improvements  followed  with  the 

development  of  the  time-motion  mode  (TM  mode)  of  display 

(27) 

and  simultaneous  recording  of  the  electrocardiogram." 

This  improvement  in  1956  by  Elder  enables  an  operator  to 
view  the  past  history  by  displaying  time  and  range  with 
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the  echo  modulating  intensity.  This  development  enabled 
the  clinician  to  quantify  structural  motion.  Although 
other  scanning  techniques  were  developed  for  pulse-echo 

ultrasonic  systems  for  medical  applications,  they  are  not 
applicable  in  cardiology  because  of  the  incompatability 
of  fast  motions  and  slow  scan  rate.  These  scanners,  such 
as  B  mode  and  compound  B  mode  scanners,  usually  produce 
transverse  cross-sections.  The  TM  mode  echocardiographic 
imaging  system  is  still  the  principal  ultrasonic  diagnostic 
tool  in  clinical  applications. 

The  application  of  B  mode  and  compound  B  mode 

(27) 

scanners  to  echocardiography  was  achieved  by  King  in  1970. 

By  using  electrocardiograms  to  synchronize  a  conventional 
B  mode  scanner,  King  was  able  to  accumulate  echos  at  the 
same  approximate  phase  of  the  cardiac  cycle  over  a  relatively 
long  period.  This  "stop  action"  technique  was  the  first 
two-dimensional  pseudo-snapshot  imaging  technique  of  the 
heart.  Various  researchers  utilized  this  technique  with 
motion  pictures  to  simulate  cross-section  versus  time 
images  of  the  heart.  This  technique  was  not  adopted  for 
general  clinical  use. 

To  image  two-dimensional  echocardiograms  in  real 
time,  scanners  were  developed  in  the  early  1970's. 

Real  time  scanners  generate  two-dimensional  images  by 
sequentially  transmitting  and  receiving  an  ultrasonic  beam 
at  different  angles,  usually  through  a  plane.  This  two 


dimensions  versus  time  technique  was  implemented  using 

( is) 

mechanical  components  by  Griffith  and  Henry.  '  By 
using  a  motor  oscillate,  the  transducer,  real  time 
sector  scans  were  provided  by  the  Griffith  and  Henry 
scanners.  This  particular  system  is  now  commercially 
available.  Electronic  means  to  sweep  the  ultrasonic 
beam  spatially  were  used  also.  By  using  multielement 
transducers  in  a  linear  array  and  activating  each 
sequentially,  Bom,  Roelandt  and  Hugenholtzv  ’  achieved  a 
real  time  planar  scan.  By  using  a  linear  phased  array, 
Thurstone^®^  was  able  to  electronically  steer  the 
ultrasonic  beam  through  a  sector  to  produce  two-dimen¬ 
sional  real  time  echocardiograms.  Such  three-dimensional 
(two  spatial  and  time)  scanners  mark  another  advance  in 
echocardiography  which  is  rapidly  being  acquired  for 
clinical  use. 

Each  of  the  newer  generations  of  echocardio- 
graphic  imaging  systems  has  advantages  over  its  pre¬ 
decessors,  and  each  system  has  its  strong  points.  For 
example,  the  earliest  A-mode  pulse  echo  system  has 
advantages  over  the  most  advanced  real-time  electronically 
steered,  phased-array  sector  scanners  in  certain 
applications.  That  is,  the  A  mode’s  high  pulse  repetition 
frequency  (PRF)  can  detect  atrial  flutter,  where  the  low 
frame  rate  of  the  sector  scanner  cannot.  In  summary,  a 
variety  of  echocardiographic  imaging  systems  exist 
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for  the  diagnosis  of  cardial  disorders  and  each  system 
has  its  own  operating  limitations.  Selection  of  one 
depends  on  the  particular  application  and  other  constraints. 

1. 2  Clinical  Applications  of  Echocardiography 

In  order  to  gain  insight  into  the  engineering 
requirements  of  echocardiography,  one  must  examine  the 
clinical  purpose  echocardiography  serves.  In  this 
subsection,  some  of  the  uses  and  proposed  uses  of  echo¬ 
cardiography  relevant  to  this  thesis  are  discussed.  The 
important  points  of  echocardiographic  diagnosis  are 
pointed  out.  These  points  are  then  translated  into 
engineering  objectives. 

Excellent  information  regarding  echocardiographic 
diagnostic  techniques  can  be  found  in  monographs  by 
Feigenbaum  and  also  in  Termini  and  Lee.^®’®®^  Numerous 
cardiac  disorders  can  be  diagnosed  by  echocardiography. 

Such  diseases  as  mitral  stenosis,  prolapsed  sintral  valve, 
aortic  stenosis,  atrial  tumors,  pericardial  effusion,  etc. 
are  routinely  detected  with  echocardiograms.  The  process 
of  cardiac  diagnosis  can  be  summarized  into  identification 
of  cardiac  structure  and  quantification  of  spatial  and 
temporal  features.  For  example,  in  the  diagnosis  of  mitral 
stenosis,  the  anterior  mitral  leaflet  is  recorded  on  a  TM- 
mode  echocardiographic  system;  then  the  velocity  of  the 
anterior  mitral  leaflet  (E-F  slope)  is  quantified  and 
compared  with  normal  values.  Here,  the  structural  dynamics 
or  temporal  features  are  exploited  for  the  diagnosis. 
Structural  dimensions  and  their  relationship  with  each 
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other  are  also  important  features  for  echocardiographic 
diagnosis.  Left  ventricular  functions  are  often  estimated 
with  TM-mode  echocardiograms.  Good  correlation  of  left 
ventricular  stroke  volume  with  cine-angiograms  can  be 
obtained  by  measuring  one-dimensional  left  ventricular 
internal  dimensional  data  at  both  systolic  and  diastolic 
instances . 

The  examples  of  assessing  the  left  ventricular 
performance  with  an  echocardiographic  imaging  system  points 
out  the  need  for  a  multidimensional  echocardiogram. 

Perhaps  more  important  than  just  assessing  the  ejection 
volume  of  the  left  ventricle  is  the  ability  of  a  multi¬ 
dimensional  imaging  system  to  quantify  local  myocardial 
hypokinesis.  This  ability  to  quantify  ischemic  dyssynergy 
is  important  in  determining  the  extent  of  coronary  artery 
disease.  Miller,  Lawther,  and  Zelis^°9^  proposed  to 
evaluate  the  area  of  ischemic  venticular  dyssynergy  by 
examining  the  regional  motion  of  the  heart  within  a  three- 
coordinate  system.  A  substitute  to  study  the  regional 
dyssynergy  without  the  use  of  a  high  data  rate  motion 
capability  can  be  achieved  by  the  comparison  of  cardiac 
wall  at  end-systolic  and  end-diastolic  instances.  Other 
established  methods  of  analyzing  ventricular  geometry 
are  biplane  ventriculographic  techniques.  "It  should  be 
noted  that  there  are  certain  aspects  of  ventriculography 
over  and  beyond  its  invasive  nature  and  the  radiation 
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exposure,  that  make  it  less  than  ideal  for  the  analysis 

of  regional  cardiac  function.”'  '  Various  researchers 

have  found  clinical  evidence  that  compound  scanning 

with  M-mode  echocardiography  has  higher  sensitivity 

in  detecting  small  regional  abnormal  motion  than  with 

(39) 

biplane  angiographic  techniques. 

In  summary,  the  development  of  echocardiographic 
techniques  has  been  on  the  increase  in  dimensions,  from 
the  one-dimensional  A-mode  to  the  three-dimensional  sector 
scanners.  An  increase  to  the  fourth  dimension  is  thought 
to  help  in  the  clinical  evaluation  of  coronary  artery 
disease  in  particular.  Although  display  techniques  are 
essentially  limited  to  three  dimensions,  it  is  important 
that  data  acquisition  be  four-dimensional;  that  is,  the 
ultrasonic  scan  should  be  accomplished  within  a  single 
heart  cycle.  This  is  especially  important  because  the 
intended  patients  might  exhibit  cardiac  arrythmia,  there¬ 
fore  rendering  data  acquisition  intervals  over  different 
heart  cycles  invalid.  It  is  toward  this  goal  of  multi¬ 
dimensional  echocardiographic  imaging  that  this  thesis 
proposes  to  reach.  In  particular,  the  specific  problem  of 
dynamic  range  compression  is  examined  and  an  algorithm  is 
proposed. 

1 . 3  Echocardiographic  Systems 

Figure  1.1  shows  the  basic  elements  in  an  echo¬ 
cardiographic  imaging  system.  The  basic  components  include 
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1)  a  transducer  to  generate  a  beam  of  ultrasonic  energy  and 
to  convert  the  reflected  acoustic  energy  into  electrical 
signal,  2)  a  transmitter  circuit  to  excite  the  transducer, 

3)  a  receiver  circuit  to  amplify  the  weak  echo  signal,  4) 
an  AGC  and/or  depth  compensation  circuit  to  compress  the 
wide  dynamic  range  of  the  echo  signal,  5)  a  bulk  storage 
or  video  processor  to  store  the  echo  signal  or  process  the 
signal  suitable  for  display,  6)  a  display  and/or  hard  copy 
device  to  present  the  image  to  the  physician,  and  7)  an 
electrocardiogram  amplifier  to  provide  synchronization 
and/or  to  provide  reference  to  the  cardiac  cycle.  Auxiliary 
subsystems  such  as  positioning  arm  with  associated  display 
circuits  or  phonocardiogram  may  be  included  in  a 
particular  display  system.  Also,  the  system  may  be  under 
computer  control.  All  echocardiographic  systems  can  be 
described  in  terms  of  the  aforementioned  configuration.  This 
artificial  characterization  eases  the  analysis  of  various 
imaging  systems.  Investigation  into  the  state  of  the  art 
of  the  various  subsystems  is  out  of  the  scope  of  this 
thesis.  Rather,  components  of  the  imaging  system  that 
affect  the  dynamic  range  of  the  echocardiographic  image  are 
discussed . 

The  basic  elements  of  the  pulse  echo  transmitter 
receiver  system  are  shown  in  Figure  1.2.  The  large  dynamic 
range  (at  least  90  dB)  of  the  echo  signal  is  the  result  of 
ultrasonic  absorption,  reflection  and  wave  divergence.  A 
grey  scale  display,  in  practice,  can  provide  about  15  dB 
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of  discernable  grey  scale  dynamic  range.  If  a  bulk 
storage  such  as  video  tape  or  digital  memory,  or  color 
display  are  used,  then  the  difference  in  dynamic  range 
may  be  less;  but  the  displayable  dynamic  range  is  much 
smaller  than  the  raw  echo  signal’s.  For  digital  memories, 
it  is  important  to  keep  the  required  dynamic  range  as  low 
as  possible  to  conserve  storage.  The  storage  constraint 
is  especially  pronounced  for  multidimensional  echocardio¬ 
graphy  because  the  amount  of  data  is  extremely  large. 

Without  loss  of  generality  for  different  echocardiographic 
imaging  systems,  the  dynamic  range  of  the  raw  echo  must  be 
compressed  at  least  two  to  three  times. 

Since  the  early  applications  of  echocardiography, 
the  problem  of  dynamic  range  compression  has  been  partially 
compensated  for  by  a  time-varying  gain.  The  basis  of  such 
an  implementation  is  that  as  the  acoustic  wave  transverses 
through  the  scanning  volume,  the  amplitude  is  attenuated. 
Therefore,  by  the  application  of  the  linearly  increasing 
dB  gain  curve  versus  range  (or  in  reality  time) ,  the 
echo  amplitude  can  be  normalized.  Operator  controls  are 
provided  so  that  the  operator  can  optimize  the  gain 
profile  for  each  patient  and  each  look  direction.  Since 
the  power  profile  for  a  ping  might  not  be  monotonically 
decreasing  at  a  single  rate,  this  technique  required  com¬ 
promises  at  higher  gains  for  certain  sections  and  saturation 
in  some  other  sections.  For  example,  echos  from  heart  valves, 
especially  at  a  high  angle  of  incidence  might  produce  low 


echos,  whereas  the  posterior  wall  produces  strong  echos. 

Compromises  must  be  made  by  the  operator  to  sacrifice 

sensitivity  within  the  atriums  or  ventricules  or  to 

sacrifice  resolution  at  the  anterior  or  posterior  wall 

(4) 

regions.  Recognizing  this  need,  Barnes  et  al.  proposed 
a  programmable  swept  gain  control.  The  basic  operating 
principle  for  such  a  dynamic  range  compression  technique  is 
through  the  use  of  voltage-controlled  gain  amplifiers  and 
wave  form  generators.  The  wave  form  generators  can  gen¬ 
erate  piecewise  linear  continuous  voltage  wave  forms  which 
can  be  used  to  control  the  gain  (in  dB)  in  the  voltage- 
controlled  amplifier.  Thus,  by  dividing  a  ping  interval 
into  multiple  segments,  the  operator  can  optimize  the  gain 
for  a  particular  function.  Also,  a  constant  rejection 
level  is  used  to  suppress  the  display  of  low  level  signals 
("grass") .  This  technique  of  multisgement  sweep  gain  has 
been  applied  to  commercial  A-mode,  TM-mode,  B-mode,  and 
spector-scanned  echocardiographic  imaging  systems.  It  was 
hailed  as  a  great  improvement  in  grey  scale  imaging. 

Although  the  manually  programmed  sweep  gain 
control  provides  great  improvement  in  dynamic  range 
compression  in  echocardiographic  imaging,  manual  adjust¬ 
ment  still  has  drawbacks.  The  quality  of  the  image 
is  still  operator  dependent,  and  nonrepeatable .  Consid¬ 
erable  time  is  required  to  set  up  the  gain  parameters 
manually.  Also,  the  gain  profile  is  constant  through¬ 
out  the  sector  scan;  thus,  for  some  look  angles,  the 
problem  of  compromising  saturation  with  low  sensitivity 
or  vice  versa  occurs.  Adaptive  or  automatic  gain  control 
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for  echocardiographic  imaging  has  been  investigated  by 

various  investigators.  DeClercq  and  Maginess  proposed 

(7) 

an  adaptive  gain  control  system  with  analog  approach. v  ' 

"The  operating  prototype  examines  peak  signals  in  eight- 
range  segments,  anticipating  required  gains  from  a  pro¬ 
gressively  updated  store.  Isolated  strong  echoes  are 
ignored  to  avoid  loss  of  lower  level  detail."  The  system 
recognizes  that  optimal  gain  occurs  over  small  spatial  or 
temporal  segments  only,  and  by  the  power  profile's  corre¬ 
lation  from  ping  to  ping,  an  updated  gain  control  is 
derived.  At  the  present,  seven  seconds  are  required  for 
the  gain  control  to  settle  for  a  constant  target.  With  a 
PRF  of  1  KHz,  that  means  seven  thousand  samples  are 
required . 

The  use  of  an  on-line  computer  for  gain  control 

( 12) 

in  ultrasonic  imaging  systems  has  been  implemented  by  Fry. 

This  system  utilizes  a  time-gated  mode.  In  other  words, 
small  segments  are  sampled  repeatedly  to  derive  its  optimal 
gain.  This  "time- gated  display"  method  requires  long  data 
acquisition  intervals  and  is  not  applicable  for  echocardio¬ 
graphy  . 

1 . 4  Scope  of  Investigation 

The  main  goal  of  the  thesis  is  to  investigate 
adaptive  dynamic  range  compression  techniques  in  echocardio¬ 
graphic  imaging  systems.  The  state  of  the  art  in  echocardio¬ 
graphy  dynamic  range  compression  is  inadequate  for  multi¬ 
dimensional  imaging.  The  manually  adjusted,  multisegmented 
gain  control  contained  in  commercially  available 
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echocardiographic  imaging  systems  is  inadequate  for  multi¬ 
dimensional  imaging  because  there  are  rapid  changes  in 
look  angles.  The  experimental  adaptive  gain  control 
proposed  by  DeClercq  and  Maginess  is  not  applicable  to 
multidimensional  imaging  because  of  the  unwieldy  amount 
of  analog  components  when  adapted  for  three-dimensional 
scanning.  Therefore,  for  multidimensional  echocardio¬ 
graphy  imaging,  an  automated,  adaptive  gain-control 
scheme  under  computer  control  is  needed.  The  use  of  a 
digital  computer  eases  the  implementation  of  the  dynamic 
range  compression  technique  and  offers  flexibility  in 
system  parameters. 

The  scope  of  this  investigation  can  be  divided 
into  three  areas.  They  are  as  follows: 

1.  An  echo-signal  model  to  provide  a  theoretical 
basis  for  the  development  of  the  dynamic 
range  compression  algorithm  and  to  provide 
simulated  data  for  the  testing  of  the 
dynamic  range  compression  algorithm. 

2.  A  general  dynamic  range  compression 
algorithm  applicable  to  any  scanning 
and  display  technique. 

3.  Examples  of  specific  implementation  technique 
of  the  dynamic  range  compression  algorithm 
and  their  performances. 

These  three  main  areas  are  covered  in  Chapters  II, 
III,  and  IV,  respectively. 


CHAPTER  II 


ECHOCARDIOGRAPHY  SIGNAL  MODEL 

The  purpose  of  this  chapter  is  to  present  an 
echocardiographic  signal  model  for  the  development  and 
test  of  the  dynamic  range  compression  algorithm.  It  is 
not  intended  to  model  every  acoustical  phenomenon  in 
the  echocardiographic  imaging  process.  Rather,  a 
statistical  approach  is  used  to  provide  a  reasonable 
signal  power  profile.  Furthermore,  the  model  is  not 
acoustically  validated  in  this  study;  as  such,  this 
model  is  based  on  published  results  from  other 
researchers.  For  the  development  and  test  of  the  dynamic 
range  compression,  an  accurate  model  of  the  echo  returns 
from  the  scanning  volume  is  not  necessary,  but  rather 
appropriate  power  profiles  are  necessary. 

Based  on  this  simplified  approach  to  modeling, 
the  process  is  illustrated  in  Figure  2.1.  The  acoustic 
wave  propagates  through  a  medium.  As  it  encounters 
reflectors,  a  portion  of  the  energy  is  reflected,  while 
the  transmitted  wave  propagates  further  into  the  medium. 
The  reflected  wave  propagates  through  the  medium  again, 
reaching  the  transducer,  and  subsequently  is  converted 
to  an  electrical  signal.  The  next  two  subsections 
examine  the  first-order  transmission  and  reflection 
processes.  Subsection  2.3  presents  a  signal  model 
with  reference  to  parameters  used  by  other  researchers. 


Figure  2.1  Signal  Model  Processes 

2 . 1  Transmission  Components 

A  single  transducer  is  used  as  the  transmitting 
and  receiving  element  in  the  pulse  mode  echocardiographic 
imaging  system.  Transducers  in  use  or  proposed  include 
unfocused  or  focused  discs  of  piezoelectric  material,  a 
linear  phased  array,  or  two-dimensional  phased  arrays. 

Often  a  mechanical  system  and  its  associated  synchronization 
electronics  are  used  with  such  a  transducer  to  form  a 
scanning  system.  In  the  case  of  phased  arrays,  this  task  is 
done  electronically.  But,  no  matter  how  the  scanning 
system  is  implemented — manually,  mechanically,  or 
electronically — the  function  of  the  transducer  is  to 
transmit  a  pulse  of  narrow  beam  acoustic  energy  into  the 
scanning  volume  and  to  receive  the  echos  produced  by 
different  anatomical  structures.  The  generation  and  pro¬ 
pagation  of  the  acoustic  wave  in  this  inhomogeneous  material 
is  a  complex  process,  but  for  practical  applications,  a 


study  of  the  first  order  effects  provides  sufficiently 
accurate  results. 

In  this  subsection,  some  of  the  factors  that 
affect  the  amplitude  of  the  acoustic  wave  are  examined. 
Here,  a  transducer  responding  to  pressure  and  operating  in 
the  linear  region  is  assumed.  Also,  the  mechanical  to 
electrical  gain  of  the  transducer  is  assumed  known  but 
ignored  in  this  discussion. 


2.1.1  Beam  Divergence 

The  spatial  distribution  of  the  acoustic  wave 
intensity  radiated  by  a  transducer  may  be  calculated  by  the 
application  ofHuygen’s  Principle.  The  intensity  pattern 
applies  to  both  transmitted  and  received  functions  of  the 
transducer  by  the  reciprocity  principle.  Since  most  echo- 
cardiographic  imaging  systems  are  designed  with  a  narrow 
beam  with  sufficiently  low  side  lobes,  and  we  seek  a  first- 
order  approximation  only,  we  will  consider  the  intensity 
pattern  along  the  direction  of  propagation  within  the 
main  lobe.  Then, for  an  arbitrary  transducer,  the  intensity 
can  be  expressed  as  a  function  of  range, 
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where  r  =  range,  along  the  direction  of  main  beam,  and 

p„T(r)  =*  temporal  maximum  pressure. 
niHx 
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Since  the  same  transducer  is  used  for  transmit  and 

receive  functions,  the  intensity  of  a  completely  reflected 

2 

wave  through  a  lossless  medium  should  be  I  (r) .  If  we  are 
interested  in  the  pressure  or  the  electrical  output  of  the 
transducer,  the  attenuation  effect  due  to  the  radiation 
from  a  transducer  is  just  I (r) . 

A  commonly  used  transducer  in  echocardiography  is 
a  one-half-inch -diameter  unfocused  transducer  operating  at 
2.25  MHz.  The  circular  disc  radiator  permits  easy  analysis 
and  is  shown  here  as  an  example  for  echo  signal  attenuation 
due  to  beam  divergence.  The  solution  to  the  vibrating 
plane  disc  is  well  covered  in  acoustic  and  ultrasonic 
imaging  literature .  Here,  the  steady  state  intensity 

distribution  is  used  for  simplicity.  Much  work  on  trans¬ 
ient  conditions  has  been  done  by  investigators  as  reported 
by  Wells. It  is  commonly  reported  that  the  solution 
converges  to  steady  state  within  three  to  six  half-cycles. 
Since  the  region  at  close  range  is  of  no  interest  (skin  and 
tissure)  to  the  cardiologist,  the  steady  state  assumption 
involves  no  degradation  to  the  model.  For  this  steady 
state  plane  disc  radiator,  the  beam  can  be  considered  to 
be  composed  of  two  regions.  Classically,  the  transition 

point  between  the  Fresnel  (near)  zone  and  the  Fraunhofer 

2 

(far)  zone  is  considered  to  be  a  /\,  where  a  is  the 
diameter  of  the  radiating  disc,  and  A.  is  the  wavelength. 

But  Fernauek's  (Ref.  70)  simulation  shows  that  the 
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transition  occurs  about  25%  closer  than  the  classical 
2 

result,  a  /A.  Nonetheless,  the  time  maximum,  spatial 
average  along  the  beam  in  the  Fresnel  zone  can  be  considered 
to  be  constant,  while  in  the  Fraunhofer  zone,  it  can  be 
considered  to  diverge.  Then, for  the  case  of  the  unfocused 
disc  transducer: 


I  (r)  = 


r  ^  ro 
r  >  r 

o 


where 


rQ  =  transition  range  between  the  Fresnel 
and  Fraunhofer  zone 

=  0.75  a 2/A  or  a2/A  (see  text). 


2.1.2  Refraction 

As  the  acoustic  beam  encounters  an  interface  where 
two  media  meets,  the  beam  is  refracted.  Section  2.2.1 
covers  the  refraction  angle  and  transmission  coefficient 
for  the  single  and  double  interface  cases.  In  summary, 
because  part  of  the  energy  is  reflected,  the  acoustic  wave 
undergoes  an  attenuation  T,  the  transmission  coefficient 
at  that  particular  range.  Also,  to  account  for  the  beam 
refracted  away  from  the  look  angle  of  the  receiver,  another 
attenuation  should  be  included.  Since  this  component  due 
to  the  refraction  angle  will  have  to  take  into  account  the 
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total  angle  as  the  beam  traverse  through  all  the  previous 
interfaces,  it  is  difficult  to  calculate  this  component. 

Since  the  velocity  of  acoustic  wave  in  biological  materials 
does  not  vary  drastically,  the  refraction  angle  is 
relatively  small  for  most  cases  and  the  exclusion  of  this 
component  is  of  no  great  consequence.  Another  term  of 
attenuation  that  occurs  at  the  boundary  is  wave  mode 
conversion.  "Mathematical  analysis  of  the  conditions 
which  arise  with  oblique  incidence  where  mode  inversion  is 
possible  is  tedious.  The  analysis  becomes  difficult  if  the 
media  are  anisotropic,  as  in  the  case  with  many  biological 
materials Therefore,  this  component  is  not  treated 
here.  The  total  effect  of  these  attenuations,  then,  is  an 
impulsive  differential  attenuation  component  occurring  at 
the  boundary . 

2.1.3  Other  Attenuation  Components 

The  reason  that  the  attenuation  components  such  as 
absorption,  scattering,  etc.,  are  grouped  together  under  one 
component  is  that  experimental  results  bear  much  more 
significance  than  theoretical  results.  Absorption  is  due  to 
the  conversion  of  ultrasonic  energy  into  heat.  When  the 
acoustic  wave  encounters  obstacles  that  are  smaller  than  the 
wavelength,  the  wave  is  said  to  be  scattered.  For  the  trans¬ 
mission  component,  we  are  interested  in  forward  scattering; 
for  reflection  component,  we  are  interested  in  back  scatter¬ 
ing.  It  is  important  to  note  that  most  biological  materials 
exhibit  frequency  dispersive  characteristics.  As  such,  in  a 
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pulse-echo  system,  the  shape  of  the  pulse  envelope  is 
changed  as  it  transverses  through  the  scanning  volume.  The 
attenuation  parameters  of  different  types  of  biological 
materials  have  been  investigated  by  various  researchers, 
and  a  summary  can  be  found  in  Wells. 

2.1.4  Total  Transmission  Effects 

As  we  have  seen  from  the  previous  subsections,  the 
total  attenuation  of  the  acoustic  beam  is  due  to  many 
processes.  In  this  subsection,  we  will  combine  these 
effects  to  model  the  transmission  process  for  our  model. 

Since  we  will  sample  the  echo  returns,  let 

A  differential  total  attenuation  in  dB's. 

Then  is  equal  to  the  sum  of  attenuation  due  to  beam 
divergence,  refraction,  absorption,  mode  conversion, 
scattering,  etc.  In  general,  one  can  approximate  the  total 
attenuation  as  a  function  of  range  for  any  scan  angle.  To 
simulate  the  actual  processes  involved,  it  is  more  practical 
to  use  measured  data  than  theoretical  results.  Furthermore, 
to  account  for  the  dispersion  effects  as  the  acoustic  pulse 
is  transmitted  through  the  scanning  volume,  one  can  modify 
the  pulse  shape  as  a  function  of  range.  Since  essentially 
no  research  data  are  available  for  dispersion,  the  property 
of  acoustic  pulse  transmission  is  ignored  in  this  simple 
model.  Also,  because  of  the  complexity  of  the  attenuation 
processes  involved,  it  is  impossible  to  account  for  all  these 


factors.  The  only  component  that  is  accountable  is  the 
beam  divergence  component.  Thus,  if  one  uses  this  to 
account  for  attenuation,  this  establishes  the  lower  bound 
on  the  attenuation  effects  as  the  acoustic  pulse  is 
traversed  through  the  scanning  volume.  This  particular 
component  is  used  in  the  dynamic  range  compression 
algorithm. 

2 . 2  Reflection  Components 

The  echo  received  at  the  transducer  is  a  summation 
of  many  phenomena  governed  by  a  complex  set  of  conditions 
as  the  acoustic  wave  propagates  through  the  scanning  volume. 

A  conceptual  understanding  of  how  echocardiography  works  is 
obtained  by  considering  the  reflection  of  a  uniform  plane 
wave  off  a  boundary  of  different  acoustic  impedance  normal 
to  the  direction  of  propagation.  But  this  simplified 
model  of  the  acoustic  process  is  inadequate.  For  example, 
as  the  acoustic  wave  transverses  through  the  scanning  volume, 
it  encounters  many  surfaces  where  the  angle  of  incidence  is 
highly  oblique.  By  Snell's  Law,  the  angle  of  reflection  is 
equal  to  the  angle  of  incidence;  then,  no  specularly 
reflected  signal  will  be  present  at  the  receiver.  But  this 
does  not  render  echocardiography  useless.  As  in  radar  and 
sonar,  backscattering  off  a  surface  is  an  important 
component.  What  the  cardiologists  are  able  to  observe  on 
existing  echocardiographic  imaging  systems  contains  a  high 
amount  of  signal  due  to  backscattering.  The  major  acoustic 
reflection  phenomena  examined  here  includes  reflection  due 
to  interfaces  and  reflection  due  to  backscattering. 
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2.2.1  Reflection  from  Interfaces 
Single  Interface 

For  a  transverse  uniform  plane  wave  incident  at 
angle  0^^  with  respect  to  normal  on  boundary  between  medium 
1  and  medium  2  (Figure  2.2),  the  angles,  reflection 
coefficient,  and  transmission  coefficient  can  be  expressed 
as  follows: 


sin  Vj 


sin 

1  N 
> 

R  = 

plr 

Pj^cosOg-pgVgCosOj^ 

pli 

P1v1cos02+P2v2cos01 

and 


p2t 

2p1v2cos01 

pli 

p Iv1cos02+P2v2cos01 

Multiple  Layer 

If  a  thin  (relative  to  pulse  envelope) 
medium  is  interposed  between  two  media,  interference  from 
the  reflection  from  the  two  boundaries  results.  Figure  2.3 
shows  such  a  configuration  for  a  normal  incidence  plane  wave. 
For  simplicity,  the  oblique  incidence  case  is  not 
discussed  here;  the  effect  of  interference  due  to  medium 
thickness  is  similar  to  the  normal  incidence  case.  The 
reflection  and  coefficients  are  as  follows: 


Figure  2.3  Reflection  From  Thin  Layer 
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and 


where 


and 


(r 


2/1 


9  9WX0 

+1>  COS^Tt(r3/X+r2/3) 


4r2/l 

9  WXo 

sin2  ( — — ) 

V* 


1/2 


r2/l 


p2V2 

P1V1 


y 


r3/l 


P3V3 

P1V1 


r2/3 


P2V2 

p3V3 


2.2.2  Scattering 

In  this  subsection,  we  consider  targets  whose 


dimension  is  smaller  than  the  wavelength.  Although  targets 
of  the  order  of  the  wavelength  are  usually  present  in  bio¬ 
logical  terms,  the  effect  of  these  scatterers  is  angle 
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dependent,  much  like  the  small  scatterer  under  discussion. 
Although  different  physical  phenomena  govern  the  two  cases, 
generality  is  not  lost  for  the  signal  model.  For  scatterers 
much  larger  than  the  wavelength,  this  is  just  a  specular 
component  discussed  under  reflection  from  boundaries. 

According  to  Rayleigh's  scattering  theory,  each 
individual  scatterer  can  be  described  by  differential 
scattering  cross-section, 


L2,/  6rKe“K 
°d  9(  X  }  a  L  K 


3Pe-3P 

2pe+p 


cos©] 


where 


\  =  wavelength, 
a  =  scatterer  radius  («A)  , 

K  =  adiabatic  compressibility  of  embedding  medium, 

Kg  =  adiabatic  compressibility  of  scatterer, 
p  =  density  of  medium 
pe  =  density  of  scatterer,  and 
9  =  angle  of  incidence. 

Note  that  is  proportional  to  the  fourth  power  of  frequency. 

(57) 

This  fact  is  confirmed  by  Shung,  et  al. 

For  forward  scattering,  one  component  of  the 
attenuation  effect,  ©=ir,  and  therefore: 
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'd  -  9^r 

9=7T 


1  /2ir\4  SrVK  3Pe'3P 


■iff)  a  L 


K  2Pe+P 


And  for  backscattering ,  9=0, 


9=0 


•i  o  4  firK  “K 
1,2ti\  _6r  e 

=  a  L-ir  + 


3pe“3p  2 
^P^P- 


Then  if  such  scatterers  are  contained  in  a  volume,  one  must 
consider  the  total  reflection  from  all  of  the  scatterers. 
Furthermore,  if  a  scatterer  moves  with  velocity  v  along  the 
acoustic  beam  direction  in  a  medium  with  velocity  of  sound 
equal  to  c,  then  this  scatterer  introduces  a  Doppler 
shift,  wi  =  2^  wQ,  to  the  carrier  frequency,  wq.  Then  to 
characterize  this  scattering  volume,  one  can  use  the 
average  scatterer  cross-section, 


ob-, f>  , 

i 

where  the  summation  is  over  all  scatterers  for  which 


and 


r .  =  r 
i 


w  ^  =  2irf  . 

For  tissue  structures,  the  velocity  of  the  scatterer  is 
relatively  slow,  then  one  can  approximate  this  as  disper¬ 
sive  in  time  only,  or 
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o(r,f)  =  <5  (f  )a(r)  . 

For  blood  cells,  where  the  individual  scatterer  velocities 
are  significant,  the  average  scatterer  cross-section  is 
both  dispersive  in  time  and  frequency. 

Furthermore  scatterers  can  be  classified  into  two 
types,  random  and  lattice.  For  randomly  located  scatterers, 
the  scattering  function  is  independent  of  the  incident  beam. 
But  for  scattering  from  a  lattice  structure,  Bragg  diffrac¬ 
tion  occurs.  Consider  adjacent  scatterers  to  be  distance 
d  apart.  By  Huygen's  principle,  the  path  length  from  the 
reflection  of  two  adjacent  scatterers  is  then  different  by 
d  sin9.  Obviously  constructive  and  destructive  interference 
can  occur  depending  on  the  difference  in  path  length.  Then 
for  scatterers  with  a  lattice  structure,  the  scattering 
cross-section  depends  on  the  angle  of  incidence. 

Shung,  et  a  1 . v  reported  that  blood  at  low 
hematocrits  (percentage  packed-cell  volume  occupied  by 
blood  cells  below  8%),  "the  red  blood  cells  behave  as 
incoherent,  randomly-distributed  scattering  centers;"  and 
"at  high  hematocrits  (greater  than  58%  or  so) ,  the  blood 

resembles  a  perfect  crystalline"  which  must  be  accounted  for. 

( 19) 

Hill  reasoned  that  "if  one  considers  the  known 
histology  of  two  human  tissues,  skeletal  muscle  and 
liver  parenchyma,  one  finds  that  they  consist  of  fibers  and 
lobules,  respectively,  whose  boundaries  are,  in  each  case, 
defined  by  connective  tissue  structures. . .  The  structures 


thus  may  constitute  arrays  which,  in  practice,  exhibit 
an  appreciable  degree  of  short  range  order...  By  assuming 
reasonable  values  for  the  parameters  describing  such 
arrays,  one  may  construct  computational  models  for  them. 

He  found  supporting  evidence  to  support  such  a  Bragg 
diffraction  theory  with  data  for  skeletal  muscle. 

In  summary,  scattering  is  a  complex  process  with 
very  few  experimental  results  where  many  parameters  are 
unknown.  If  we  consider  the  effect  of  the  scatterers  at 
one  point  and  at  a  single  angle,  we  can  assume  that  it  is 
a  Rayleigh  process.  But  the  parameter  is  a  difficult  to 
predict  function  of  many  factors.  Also,  frequency 
dispersion  introduces  a  time-varying  filter  property  to 
the  pulse  envelope.  Dispersion  is  unimportant  to  our 
signal  model  because  the  effect  is  greatest  for  blood  cells, 
which  are  of  little  interest  in  an  echocardiographic  image. 

2 . 3  Signal  Model 

The  preceding  sections  discussed  some  of  the 
mechanisms  which  affect  the  echo  returns  in  echocardiography 
Since  echocardiographic  imaging  requires  information  through 
out  the  scanning  volume  to  generate  an  image,  we  are  con¬ 
cerned  with  phenomena  taking  place  throughout  the 
scanning  volume.  This  particular  aspect  of  echocardio¬ 
graphic  imaging  is  dissimilar  to  radar  or  sonar  analysis. 

In  the  latter,  the  volume  of  interest  spans  only  a  small 
part  of  the  scanning  volume.  Even  for  the  case  of  multiple 
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targets,  for  example,  aircraft  in  formation,  the  problem 
is  mathematically  tractable  because  of  the  target 
similarity  and  of  the  propagation  medium.  In  echocardio¬ 
graphy,  the  process  is  much  more  complicated  because  of 
the  medium  and  of  the  continuum  of  highly  dissimilar 
targets.  Therefore,  a  simplified  signal  model  is  presented 
here  for  simulating  the  power  level  profile  rather  than  the 
acoustic  processes  involved. 

2.3.1  Reflection  from  a  Small  Volume 

From  the  preceding  subsections,  we  see  that  the 
received  signals  are  composed  of  returns  from  distinct 
fixed  targets,  such  as  boundaries  of  acoustic  impedance  mis¬ 
match  and  from  point  targets,  such  as  random  or  structured 
scatterers.  Then,  without  lost  of  generality,  we  can  model 
the  echo  as  a  sum  of  a  specular  and  a  scatterer  component. 

(If  only  the  specular  component  is  present,  one  can  set  the 
scatterer  component  to  zero,  or  vice  versa.)  In  Section 
2.2.1,  we  saw  that  the  reflection  from  the  interface  between 
two  media  is  highly  angular  dependent.  There  are  probably 
more  oblique  angles  of  incidence  than  occurrences  of  normal 
incidence.  (Consider  a  sphere  or  ellipsoid,  for  instance; 
from  any  scanning  angle,  there  are  more  highly  oblique  angles 
of  incidence  than  those  close  to  normal  incidence.)  Also, 
for  irregular  interface  surfaces,  such  as  the  endocardium, 
there  are  different  angles  to  incidence  within  the  beam 
width.  Therefore,  there  is  probably  a  preponderance  of 
scatterer  returns  rather  than  specular  returns. 
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For  simplicity  in  this  discussion,  let  us  ignore 
the  time  delay  involved.  That  is  if  pulse  p(t)  is  trans¬ 
mitted  at  t=0,  the  echo  from  a  target  distance  x  away  will 

be  received  as  p(t-2x/v„) ,  where  v„  is  the  harmonic  velocity 

rX.  n  n 

and  is  equal  to  Aj  ^yydy j  1.  Then,  ignoring  the  delay,  a 
general  received  signal  from  a  single  target  and/or  small 
scatterer  volume, 

r(t)  =  ap(t)  cost w  t  +  0  (t)  +  d] 

c  p 

+  Bp(t)  coslw  t  +  0  (t)  +  0], 

C  P 

0  <  t  <  T  . 

where : 

r(t)  =  echo  signal  received  at  the  transducer, 

p(t)  =  transmitted  pulse  envelope, 

0p(t)  =  phase  of  transmitted  pulse, 

w  =  carrier  frequency, 
c 

a  =  amplitude  of  the  specular  component, 

6  =  phase  of  the  specular  component, 

3  =  amplitude  of  the  scatterer  component,  and 

0  =  phase  of  the  scatterer  component. 

Since  echocardiography  pulse  echo  systems  are  not  phase 
modulated,  0p(t)=O.  In  this  case,  a  and  6  are  assumed  to 
be  known.  Without  loss  of  generality,  we  can  let  6=0. 
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Decomposing  r(t)  into  the  quadrature  terms, 

r(t)  -  al(t)Lp(t)coswct]  +  ag (t) [ p(t) sinwct] , 

where  a^  and  a2  are  independent  Gaussian  random  variables 
with 

ECa^  =  a  , 

E(a2)  =  0  , 

and 

Var(a1)  =  Var(a2)  =  a2 

Then,  at  the  output  of  the  detector,  the  envelope  is 
s(t)  =  Ja2  +  a2  p(t) 

The  probability  density  function  for  s(t)  is  Rician, 

(s2+q) 

P,(»  -  -Tf  e  2o*  V2!’  u<s>’ 

a  a 

where  IQ(x)  -  Bessel  function  of  the  first  kind. 
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Actually,  the  above  describes  three  basic  situations: 
reflection  from  a  boundary,  reflection  from  a  heterogeneous 
tissue  volume,  and  reflection  from  a  blood  volume.  Con¬ 
ceptually,  we  desired  to  identify  each  of  the  three 
situations  and  to  estimate  the  parameters  involved.  The 
parameters  for  the  different  situations  are  from  different 
underlying  populations.  For  the  reflection  from  a  boundary, 

reflection  coefficient  (normalized  transmit 
pulse) 

backscatter  component  of  the  boundary  surface. 

For  the  reflection  from  heterogeneous  tissue  volume, 


a 


and 


and 


a  =  0 

o 

o  =  backscatter  component  of  the  tissue  volume. 


For  the  reflection  from  blood  volume, 


and 


a  =  0 

2 

o  =  backscatterer  component  from  blood  cells. 


Or,  from  the  detection  viewpoint,  there  are  two  hypothesis: 
Hp  specular  component  present;  and  HQ,  only  scatterer 
component  present.  The  HQ  hypothesis  can  be  further  broken 


35 


2 

down  for  high  and  low  o  ;  the  latter  with  a  high  probability 
of  being  a  volume  of  blood  only. 

2.3.2  Complete  Signal  Model 

If  we  are  to  consider  the  total  echo  return  as  the 

acoustic  wave  traverses  through  the  scanning  volume,  we  have 

to  modify  the  aforementioned  small  volume  model.  This  can 

be  accomplished  by  first  assuming  the  pulse  envelope  to  be 

sufficiently  short  in  duration  compared  with  the  sampling 

time.  Then,  the  sampled  envelope  signal  for  a  target  volume 

at  distance  x .  is 
J 

s(i)  -/■;«>  +  a^(j)  6(i-j)  , 

where 


s(i)  =  echo  detector  output  at  the  ith  sampling 
time, 

a0  =  Quadrature  components  due  to  reflection 

’  from  target  volume  at  x.,  and  processes 

statistically  properties  as  in  Sub¬ 
section  2.3.1, 


<5 ( i— j )  =  Kronecker  delta 


sample  delay  corresponding  to  distance  Xj 

(-  -v^.-  ,  where  At  =  sampling  interval)  . 

AtVH 
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Suppose  we  introduce  a  pulse  function  p(t)/0  for  t  At, 
and  with  pulse  samples,  p(K)=p(tK),  for  k=0,l,...,K,  then 
the  quadrature  components  will  have  to  be  modified  to 

i 

^p(i-k)ao  l(k)  , 

k=i-K 

or,  since  p(k)=0  for  k>K, 

i 

y^p(i-K)ao>1(K)  , 

k=irK 

Now,  if  we  introduce  the  total  attenuation,  we  have  to 
introduce  a  multiplier,  m(k) ,  where 

k 

m(k)  =7 1  a0  » 

^=0  x 

and 


a D  =  total  attenuation  per  sample  distance. 


Then  the  detector  output  becomes 


s  (i) 


Figure  2.4  shows  a  block  diagram  of  the  signal  model. 
Note  that  the  process  is  still  an  envelope  of  the  quad¬ 
rature  terms, 
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But  now  with  the  introduction  of  the  pulse  envelope  and 
total  attenuation  terms,  the  quadrature  terms  become 

i 

bQ  ( i)  =  m(k)aQ  1(k)p(i-k)  . 
k=i-K 

Since  a  .  is  a  normal  process,  then  b  ,  remains  a  normal 
o,l  o,l 

process,  but  means  and  variances  change  to  the  following: 


and 


i 

Eih0>^  =  y^m(k)p(i-k)E^a0t  1(k)| 
k=i-K 

i 

°2lb o,J  =  °2[m(k)p(i-k)ao  l(k)J 

k=i-K 


i  i 

+  2  Cor^m(k)  p(i-k)aQ  j^k), 

k=i-K  A=i-K 
k  tSL 

n(i)p(H)a0  AH)}  • 


Note  that  even  if  a  , (k)  are  uncorrelated  in  range,  with 

O  ,  I 

the  introduction  of  m,  p,  the  modified  quadrature  terms 
will  be  correlated  in  range.  For  each  sample,  though,  it 
still  has  a  Rician  distribution.  The  variance  for  such  a 


Figure  2.4  Echocardiography  Signal  Model 
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sample  is  now  increased  due  to  the  correlated  nonstationary 

/“ o  T~ 

ao+ai  , 

we  would  desire  to  make  the  process  stationary.  Choosing 
p(t)  short  is  one  step  toward  avoiding  correlation 
induced  by  the  ensonif ication .  This  limits  the  available 
power,  especially  in  medical  applications  where  peak 
power  must  be  constrained.  Pulse  compression  techniques 
such  as  those  utilized  in  radar  might  be  considered  here, 
but  the  complexity  and  expense  preclude  the  feasibility 
of  such  systems  at  that  time.  The  other  approach  to  making 
the  processes  stationary  then  is  to  compensate  for  m. 
Anyway,  it  is  necessary  to  compress  the  dynamic  range  of 
the  signal  into  the  limited  dynamic  range  of  a  real  pro¬ 
cessor.  The  effect  of  eliminating  m  then  is  equivalent  to 
using  an  even  illumination  function,  as  if  the  transducer 
is  directly  in  front  of  each  target.  It  seems  reasonable 
to  estimate  m  and  divide  by  the  estimate.  Since  m  is  a 
complicated  process  whose  effects  show  up  within  the  signal 
information  itself,  a  direct  estimation  is  infeasible. 

The  dyanmic  range  compression  algorithm  to  be  presented 
contains  a  multiplier,  which  is  proportional  to  the 
received  power.  Over  small  volume,  this  represents  a  rough 

estimate  of  m  and  a  .  Therefore,  it  tends  to  normalize 

ao,l 

the  signal  power  level  throughout  the  scan  volume. 

2.3.3  Example 

In  Subsection  2.3.2,  a  model  of  the  complete  echo 
return  along  a  scan  direction  was  presented.  In  this 
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subsection,  examples  of  how  the  parameters  for  the  model 
can  be  obtained  are  given.  As  mentioned  earlier,  the 
acoustic  reflection  and  propagation  processes  are  immensely 
complex,  the  method  only  serves  to  estimate  power  profiles 
rather  than  a  complete  acoustic  model.  Therefore, 
simplifications  were  made  where  the  power  profile  is  not 
severely  affected. 

Figure  2.5  shows  a  typical  single  scan  in  echo¬ 
cardiography.  The  region  of  interest  spans  from  the 
anterior  epicardium  to  the  posterier  chest  wall. 

In  the  chest  wall,  the  ultrasound  is  mainly 
scattered  by  the  muscular  tissue.  Actually,  a  complicated 
phenomenon  occurs  in  this  near  zone.  Reverberation  is 
another  possibility.  Also,  the  sternum  and  ribs  tend  to 
contribute  dislocation  from  the  ideal  case  also.  Very 
little  experimental  work  has  been  done  in  this  area.  Since 
this  region  is  of  no  interest,  the  simple  model  of  just 
scattering  is  sufficient. 

As  the  ultrasonic  beam  traverses  through  the  scanning 
volume,  it  encounters  tissues  and  fluids  of  different 
acoustical  impedance.  The  reflection  coefficient  and 
transmission  can  be  calculated  at  these  interfaces  from 
their  characteristic  impedances  and  angle  of  incidence. 

Since  the  transducer  serves  as  the  transmitter  and  the 
receiver,  only  the  reflected  component  in  the  direction  of 
incidence  should  be  considered.  This  can  be  estimated  by 
multiplying  by  the  transducer's  directivity  function;  thus, 
the  specular  component  can  be  estimated.  At  the  same  time, 


Figure  2.5  Example  of  Signal  Model  Parameter 
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the  transmission  coefficient  on  the  return  trip  can  be 

calculated.  The  two  transmission  coefficients  can  then 

be  combined  with  the  attenuation  component  for  that 

medium  to  form  the  total  attenuation  parameter.  The 

attenuation  and  scatterer  component  parameters  at  the 

different  ranges  can  be  set  according  to  the  value  for  the 

medium.  The  scatterer  component  at  the  interfaces  has  to 

be  changed  to  account  for  scattering  off  the  surface. 

Since  an  overall  power  profile  is  desired,  it  is  not 

necessary  to  adjust  the  scatterer  component  according  to 

the  angle  of  incidence.  An  example  of  parameters  is 

illustrated  in  Figure  2.5. 

The  velocity,  absorption,  attenuation  parameters 

and  characteristic  impedances  have  been  investigated  by 

f  68 ) 

many  researchers  and  is  well  summarized  in  Wells.  ' 

For  echocardiography,  a  transducer  frequency  of  2.25  MHz 

is  normally  used.  For  this  frequency,  the  parameters  are 

/ 

summarized  in  Table  2.1. 


Table  2.1  Ranges  of  Parameters  in  Biological  Mater i 
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CHAPTER  III 


DYNAMIC  RANGE  COMPRESSION  TECHNIQUES 

The  object  of  the  section  is  to  develop  a  general 
dynamic  range  compression  algorithm  applicable  to  any  of 
the  scanning  and  display  techniques.  Although  the  hard¬ 
ware  implications  of  various  echocardiographic  imaging 
systems  are  drastically  different  from  each  other,  the 
basic  pulse  echo  technique  is  common  to  all.  The 
receiver  turns  on  for  a  period  of  time  before  the  trans¬ 
mitter  can  be  activated  again.  This  ping-by-ping  or 
line-by-line  approach  is  common  to  all  of  the  scanning 
display  methods.  Only  the  transducer  orientation  and 
display  format  are  different.  In  this  section,  the 
organization  of  echo  signal  for  various  systems  is  defined. 
Based  on  the  statistical  echo  signal  model  suggested  in 
Chapter  II,  a  general  dynamic  range  compression  algorithm 
for  echocardiographic  imaging  systems  is  proposed.  Basic 
implementation  and  feasibility  of  the  general  algorithm 
are  examined. 

3 . 1  Introduction 

Although  application  of  dynamic  range  compression 
techniques  to  echocardiographic  imaging  has  been  limited  as 
shown  in  Section  1.3,  a  wealth  of  information  on  automatic 
gain  control  and  dynamic  range  compression  in  the  radio, 
radar,  sonar,  and  data  transmission  area  exists.  Also,  the 
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process  of  compressing  the  dynamic  range  of  an  echocardio¬ 
gram  is  an  image  processing  technique.  Therefore,  relevant 
investigation  3n  this  area  exists.  There  are  special 
features  and  objectives  where  these  well-explored  off-line 
computer  techniques  cannot  be  applied  to  echocardiographic 
image  processing  directly,  but  these  results  have  a 
significant  impact  on  the  algorithm  proposed  here. 

3.1.1  Automatic  Gain  Control  in  Communication 


Automatic  gain  controls  (AGC)  are  used  extensively 
in  receiver  circuits  to  minimize  gross  variation  on  signal 
power  levels.  The  most  common  purpose  of  an  AGC  system  is 
to  suppress  low-frequency  variation  in  signal  strength 
without  disturbing  the  higher  frequency  signal  components. 
Numerous  AGC  configurations  have  been  designed  to  accomplish 
this  goal.  For  simplicity,  AGC  systems  can  be  considered  to 
be  one  of  three  different  types:  feedback,  feedforward  and 
instantaneous  types.  Various  configurations  of  AGC  circuits 
were  adopted  for  use  in  echocardiography.  Since  these 
circuits  are  generally  used  to  suppress  low  frequency  power 
level  variation,  they  essentially  process  in  one  dimension 
only.  Therefore,  these  AGC  systems  must  be  modified  to 
some  extent  to  accommodate  multidimensional  imaging. 

A  commonly  used  AGC  system  in  echocardiography  is 
the  implementation  of  instantaneous  automatic  gain  control 
(IAGC)  with  logarithmic  amplifiers.  These  IAGC’s  might  or 
might  not  be  used  with  sweep  gains.  The  difference  is  that 
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their  compression  ratios  are  different.  The  dominant 
appeal  of  such  logarithmic  amplifiers  in  the  echo  receiver 
chain  is  that  a  human  responds  to  logarithmic  levels 
visually.  This  definitely  serves  as  a  matching  transformer 
between  the  electrical  signal  with  the  human  visual  response 
system.  But  if  the  signal  is  stored  in  some  kind  of  bulk 
storage,  for  example,  a  digital  memory,  then  the  question  of 
matching  to  the  recording  medium  arises.  The  logarithmic 
compression  to  match  the  human  response  can  always  be 
incorporated  into  the  display  system,  and  the  signal  should 
then  be  matched  to  the  recording  medium.  In  terms  of 
quantization,  if  digital  memory  is  used,  logarithmic 
compression  has  tremendous  appeal.  Signals  quantized 
uniformily  after  logarithmic  amplification  provide  constant 
percentage  error.  This  uniform  percentage  error  property 
is  valid  for  any  signal  distribution.  But  for  human 
distribution,  logarithmic  quantization  is  not  optimal  in 
the  least  square  error  sense.  Also,  for  echocardiographic 
imaging  applications,  the  logarithmic  compression  ratio 
must  be  set  over  the  whole  scanning  volume;  thus  for  local 
neighborhoods  with  fairly  constant  power  profiles,  the 
compression  ratio  cannot  be  set  optimally  in  order  to 
satisfy  the  global  requirement.  Therefore,  IAGC's  are 
inadequate  for  echocardiography  dynamic  range  compression. 

Radio  type  feedback  automatic  gain  controls  are 

(34) 

used  in  some  echocardiography  instruments.  Figure  3.1 

(7) 

illustrates  such  an  AGC  system.  DeClercqv 


summarized 


Figure  3.1  Radio  Type  Feedback  Gain  Control 
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the  gain  control's  inadequacy  as  follows: 


"[However,]  this  method  is  basically  inappropriate 
to  the  type  of  signal  in  ultrasonic  pulse  echo 
systems.  Conventional  radio  and  television 
signals  are  modulated  onto  a  continuous  carrier. 
Detection  provides  an  average  output  readily 
extracted  with  the  simple  filters  indicated. 
Ultrasound  signals  are  discontinuous  with 
periods  of  zero  return  over  which  the  best  that 
can  be  achieved  is  gain  increase  until  the 
noise  output  rises.  When  the  next  distinct  echo 
occurs,  receiver  gain  can  start  to  adjust.  If 
adjustment  is  fast  enough  to  follow  the  fastest 
pulses  the  system  becomes  unstable.  Longer 
response  times  permit  transitory  overloads  and 
the  final  gain  value  reached  depends  on  the 
duration  as  well  as  the  amplitude  of  the  echo." 

DeClercq  further  suggested  that  a  dual  channel  feed¬ 
forward  configuration  such  as  Figure  3.2  can  overcome,  to 
some  extent,  the  limitations  in  the  feedback  configuration. 
"This  arrangement  still  suffers  from  the  basic  limitation 

that  it  acts  with  only  local  knowledge  of  the  signal 
(7) 

statistic."'  '  Based  on  this  configuration,  DeClercq  and 
Maginess  developed  the  adaptive  gain  control  for  dynamic 
ultrasound  imaging. 

In  summary,  radio  type  automatic  gain  controls  are 
inadequate  for  multidimensional  imaging  limitations  in  taking 
advantage  of  the  spatial  and  the  temporal  correlations  which 
exists  in  multidimensional  images. 


3.1.2  Adaptive  Gain  Control  for  Dynamic 
Ultrasound  Imaging 

( 7 ) 

DeClercq  and  Maginess'  disclosed  an  adaptive  gain 
control  for  dynamic  ultrasound  imaging.  A  schematic  for 
such  a  system  is  shown  in  Figure  3.3.  "The  operating 


Figure  3.3  Adaptive  Gain  Control  System  (DeClercq  and  Maginess 


prototype  examines  peak  signals  in  eight-range  segments, 

anticipating  required  gains  from  a  progressively  updated 

store.  Isolated  strong  echoes  are  ignored  to  avoid  loss 

(7) 

of  lower  level  detail."  This  feedforward  adaptive  gain 
control  utilizing  temporal  (if  manually  scanned,  spatial 
also)  correlation  of  a  power  profile  has  been  tried  on  an 
M-mode  echo  cardiographic  imaging  system.  It  was  shown 
that  a  satisfactory  gain  profile  was  reached  in  seven 
seconds  for  a  constant  target.  Apparently,  this  is  satis¬ 
factory  if  M-mode  images  are  desired. 

This  analog  technique  is  an  automated  adaptation  of 
sweep  gain  control,  which  has  clinically  shown  its  grey  scale 
capability.  The  long  settling  time  or  learning  period 
precludes  real  time  or  fast  scanning  at  wide  look  angles. 

Its  adaptability  shows  its  main  advantage  in  adjusting  to 
different  patients.  Due  to  its  long  settling  time,  this 
technique  is  inadequate  for  multidimensional  echocardio- 
graphic  imaging. 

Other  practical  considerations  hamper  this 
technique  in  multidimensional  imaging.  For  example,  for 
sector  scanning,  if  different  gain  profiles  are  used  for 
small  scanning  angles,  then  banks  of  capacitors  are 
required  for  memory.  An  extension  to  three  or  four  dimen¬ 
sions  will  definitely  be  uneconomical  to  implement.  The 
large  number  of  "memory"  capacitors  required  is  further 
increased  if  the  number  of  segments  is  increased  to  16, 
which  is  commercially  available  in  manually  adjusted  sweep 
gain  echocardiographic  systems. 


The  problem  of  the  large  amount  of  analog  memory 
required  can  be  solved  by  the  use  of  digital  memory  and 
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the  use  of  digital  computer.  This  increases  the 
flexibility  of  the  control  scheme.  The  use  of  a  digital 
computer  permits  statistical  techniques  rather  than  a 
"heuristic"  implementation  in  the  analog  system.  Also,  the 
traditional  constant  rejection  level  can  be  implemented  in 
the  same  fashion.  These  value-added  features  form  the 
basis  of  the  dynamic  range  compression  algorithm  presented 
in  this  thesis. 


3.1.3  Statistical  Techniques  for  Gain 
Control 

A  gain  control  problem  similar  to  that  in  echo- 

cardiographic  imaging  systems  exists  in  communication  and 

radar  receivers.  The  statistical  aspects  of  gain  control 

(47 

in  such  applications  is  well  explored  in  the  literature. 

51) 

The  most  applicable  situation  can  be  found  in  a  paper 
(49) 

by  Rappaport.  The  problem  of  optimum  AGC  in  an  inco¬ 

herent  receiver  is  explored. 


"A  Bayesian  learning  algorithm  for 
statistically  optimum  adaptive  gain  control 
(AGC)  is  derived  in  which  the  receiver 
gain  is  adjusted  during  the  learning 
process  so  that  the  average  cost  of 
excluding  the  received  signal  from  the 
receiver's  dynamic  range  is  minimized. 
Application  of  the  approach  to  the 
frequently  occurring  situation  of  Raleigh- 
distributed  signals  of  unknown  strength 
yields  a  recursive  algorithm  that  can  be 
easily  implemented. " (’“) 
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Although  these  techniques  form  a  good  statistical 
basis  for  a  statistically  optimum  gain  control  to  minimize 
exclusive  cost,  it  lacks  applicability  to  echocardiography 
in  certain  aspects.  In  echocardiographic  images,  the 
"clutter"  consists  of  scatterer  returns  from  blood,  and  any 

a 

signal  above  the  clutter  level  is  of  interest  to  the 
operator.  Therefore,  the  lower  limit  should  be  set  by  the 
clutter  level  rather  than  by  considering  the  exclusion  cost 
as  suggested  in  the  algorithm.  An  improvement  which 
eliminates  the  learning  nature  of  the  algorithm  can  be 
accomplished  by  using  a  separate  power  level  estimation 
channel.  These  statistical  techniques  must  be  modified  for 
echocardiographic  imaging  applications. 

3.1.4  Image  Processing 

The  dynamic  range  compression  problem  in  multi¬ 
dimensional  echocardiography  can  be  considered  as  an  image 
processing  problem.  Consider  any  two  dimensional-plane  of 
the  multidimensional  echocardiographic  image  (for  example, 
the  M-mode  or  the  B-mode  images) .  The  object  is  to  improve 
the  grey  scale  assignment  of  such  an  image.  This  is  also 
equivalent  to  produce  an  even  illumination  function.  This 
problem  and  other  associated  image  processing  techniques 
have  been  investigated  generally  and  specifically  for  echo¬ 
cardiography. 

Various  investigators  have  implemented  echocardio¬ 
graphic  image  processing  techniques  incorporating  a  digital 

(3) 

computer.  Ballard,  Waag  and  Angel  v  use  two-dimensional 
Fourier  Transform  techniques  to  enhance  echocardiographic 
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images  by  high  and  low  pass  filtering.  Waag  and  Gramiak^66^ 
use  digital  linear  filtering  to  distinguish  and  separate 
signal,  noise,  and  clutter  components  in  an  augmentation  of 
M-mode  scanning.  Other  image  processing  techniques  besides 
amplitude  manipulation  have  been  investigated  also. 

McSherry  and  Keller,  ’  for  example,  incorporate  by  means 
of  a  digital  computer  an  inverse  filter  for  deconvolving 
the  effects  of  the  ultrasonic  beam  width.  Digital  inter¬ 
polation  in  both  time  and  distance  simultaneously  have  been 
used  in  an  attempt  to  eliminate  artifacts,  to  diminish  the 

effects  of  drop-out  with  unfavorable  angles  of  reflection, 

(35  63-67) 

and  to  facilitate  data  expansion.  ’  '  Although 

these  techniques  show  promise,  all  of  the  data  have  been 
digitized  from  commercial  echocardiographic  units  or 
through  video-tape  recorders.  The  dynamic  range  of  their 
data  has  been  compressed  by  sweep  gain  and  rejection  cir¬ 
cuits.  The  digitizers  quantize  the  data  to  eight  bits, 
equivalent  to  48  dB.  Therefore,  even  though  these 
techniques  show  promise  in  compressing  the  48  dB  data  to 
the  display  level  of  about  24  dB,they  are  not  applicable 
in  real  time  dynamic  range  compression  where  the  raw  echo¬ 
cardiographic  signal  must  be  compressed  to  display  dynamic 
range . 

Techniques  for  general  image  processing  have 

implications  in  the  echocardiography  dynamic  range  com- 

(1  17—18) 

pression  problem.  "Andrews,  Hall  and  othersv  ’  ’  have 

produced  enhanced  imagery  by  a  histogram  equalization 
process  for  which  the  histogram  of  the  enhanced  image  is 
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forced  to  be  uniform. To  minimize  the  effects  of 

nonuniform  illumination,  Oppenheim,  Schaffer  and 

(44) 

Stockhamv  suggest  the  use  of  homomorphic  filter  with 
a  high  pass  function.  These  two  techniques  have 
applications  in  echocardiographic  image  processing,  but 
they  are  post-processing  techniques.  Again,  a  wide 
dynamic  range  image  data  must  be  quantized  for  computer 
processing,  thus  the  methods  cannot  be  utilized  for  echo¬ 
cardiographic  dynamic  range  compression. 

3 . 2  Development  of  Algorithm 

It  is  desired  to  develop  a  general  dynamic  range 
compression  algorithm  applicable  to  multidimensional  echo¬ 
cardiographic  imaging.  The  algorithm  should  be  implemented 
in  or  near  real  time  under  computer  control.  It  should  be 
capable  of  compressing  the  echo  level  dynamic  range  (~90  dB) 
to  a  display  level  dynamic  range  («vl5  dB)  or  to  a  storage 
level  dynamic  range  (example,  48  dB) .  The  latter  permits 
flexibility  for  post-processing;  but  the  same  algorithm 
can  be  applied  again  to  reduce  the  dynamic  range  to  display 
level.  Without  loss  of  generality,  a  grey  scale  display  is 
considered  because  color  assignment  is  arbitrary  and  can 
be  considered  a  post-processing  technique. 

3.2.1  Objectives 

The  objectives  of  producing  a  gr.ey  scale  echocardio¬ 
graphic  image  can  be  summarized  as  follows: 
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MINIMIZE: 

A.  Effects  of  attenuation,  absorption, 
divergence,  etc.  In  other  words, 
provide  a  constant  illumination  function. 

B.  Clutter.  Extraneous  returns  such  as  back- 
scattering  from  a  blood  volume  should  be 
eliminated  from  the  display. 

C.  Artifacts  due  to  the  AGC  algorithm.  The 
algorithm  should  not  generate  false 
structures . 

D.  Saturation.  Although  saturation  at  isolated 
locations  does  not  degrade  the  image  full 
brightness  over  a  local  neighborhood  it  tends 
to  degrade  resolution. 

MAXIMIZE: 

E.  Detection  of  structures  at  high  angles  of 
incidence.  Cardiac  structures  such  as 
values  leaflets  are  of  great  interest,  but 
they  usually  are  at  high  incidence  angles  to 
the  ultrasonic  beam.  Elimination  of  this 
"wash-out”  effect  is  highly  desirable. 

F.  Use  of  dynamic  range  of  display.  Full 
dynamic  range  of  the  display  unit  should 
be  utilized  over  a  local  neighborhood. 

3.2.2  Definitions 

Although  there  are  many  types  of  scanning  and  display 
methods  in  echocardiographic  imaging,  dynamic  range  compression 
or  signal  processing  problems  are  common  to  all.  The 
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difference  between  these  different  types  lies  in  the 
geometry  and  time-frame  organization.  It  is  intended 
in  this  subsection  to  unify  the  difference  scanning  and 
display  methods  so  that  the  general  dynamic  range  compression 
algorithm  is  common  to  all  but  the  actual  implementation 
method  varies. 

The  objective  of  echocardiography  is  to  visualize 
a  solid  and  its  motions.  Thus,  echocardiographic  imaging 
systems  are  organized. so  that  one  or  more  of  the  four 
cardinal  dimensions,  the  three  geometric  dimensions  plus 
time,  can  be  viewed  alone  with  the  echo  amplitude.  For 
example,  for  M-mode  devices,  range  and  time  are  the  indepen¬ 
dent  variables  and  the  echo  amplitudes  is  a  function  of 
these  two  variables.  Perhaps  at  this  stage,  the  variables 
range  and  time  should  be  clarified.  Range  is  actually  a 
time  function,  since  it  is  proportional  to  the  product  of 
sound, velocity  and  time.  For  echocardiography,  the 
maximum  range  of  interest  is  usually  about  fifteen  centi¬ 
meters  with  the  corresponding  maximum  time  equal  to  about 
two-hundred  microseconds.  Since  the  cardiac  structure 
moves  with  a  much  slower  period,  this  is  essentially  a  snap¬ 
shot  of  a  line  of  tissues.  The  time  variable  as  defined 
here  refers  to  the  next  group  of  transmit-and-receive 
cycles,  or  the  neat  frame.  With  these  definitions  for  the 
geometric  and  temporal  dependencies,  one  can  combine  all 
the  different  types  of  scanning  and  display  methods  and 
define  the  received  echoes  or  other  variables  as  in  Table 
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Table  3.1  Independent  Variables  for  Different 
Scanning  Modes 


Mode 

Continuous 

Discrete 

M  mode 

Sector  Scan  Snap  Shot 

Real  Time  Sector  Scan 

3D  Snap  Shot 

3D  Real  Time 

(r.  \£) 

(r,  e) 

(r,  e,tf ) 

(r.e,<t>) 

(r,  e,  <D,ty  ) 

(rj(t^)  or  ( i,  £  ) 

(rj(e.)  or  ( i,  j) 

(rj,  ej,t  jp)  or  (i,j,f) 

(rj,0j.  <»k)  or  (i,  j,  k) 

(ri'ej’  ‘V* l]  or  (i’  ^  k’^  1 

3D  Real  Time 


The  difference  for  different  scanning  and  display 
methods  is  illustrated  in  Figure  3.4.  The  data  acquisition 
intervals  are  illustrated  in  the  figure.  Times  for  display 
can  be  concurrent  with  the  data  acquisition  interval  or 
bulk  storage  can  be  used  and  then  the  image  displayed  at 
the  end  of  each  frame.  Without  loss  of  generality,  the 
independent  variable  defined  on  Table  3.1  apply  to  each  of 
these  scanning  methods  and  (•)  denotes  these  independent 
variables  depending  on  the  method  used. 

The  echo  signal  and  other  variables  are  defined  as 

follows : 


s(-)  =  received  echo  amplitude, 

v( • )  =  dynamic  range  compressed  video, 

vq(')fq  =  1 , 2 , . . . , Q]  =  quantized  v(-) 

or  discrete  grey  or  color  level  assignments, 
g(’)  =  gain  level, 
b(-)  =  bias  level, 

I  =  number  of  echo  amplitude  samples  per  ping, 

J  =  number  of  pings  per  sector, 

K  =  number  of  sectors  per  frame, 

L  =  number  of  frames, 

M  =  IxJxKxL  =  total  number  of  echo  amplitudes,  and 
Sp  =  set  of  all  possible  subset  of  data  points 

=  [spv  Sp2 . ~prJ  > 


where 


Figure  3.4  Timing  Diagram  for  Different  Modes  of  Echocardio 
graphic  Imaging 
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spm  =  all  possible  subsets  of  data  points  with 
m  points  in  subset, 
p  =  index  for  subset  =  1,2,...,P,  and 
m  -  number  of  points  in  subset  =■  1,2 . M  . 

Therefore,  if  there  are  M  data  points, 

M 


3.2.3  Approach 

The  objective  of  determining  volumes  of  similar 
biological  material  and  ascertaining  their  attributes  has 
a  direct  theoretical  solution.  If  one  takes  all  combinations 
of  the  received  echo  amplitude  samples  and  uses  them  to  test 
the  hypothesis  of  samples  from  a  certain  population,  then 
one  solves  the  problem  of  clustering  samples  from  the  same 
biological  material  and  estimating  its  parameters.  This 
method  will  lead  directly  to  a  likelihood  ratio  test: 

H 

A(V  2  *  • 

where 

H  =  hypothesis  that  these  group  of  echoes  belong 
to  similar  material, 

H  =»  negative  of  H,  and 
rj  -  threshold  value. 
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m 

Since  there  are  P  =  (M)  combinations  of  s  and 

l—u  m  p 

m-1 

M  -  IxJxKxL,  this  approach  is  impractical;  P  is  an  astro¬ 
nomical  number.  Even  if  we  limit  the  above  combinations 
to  sets  that  make  geometrical  sense  (neighborhood  echo 
amplitudes  have  a  high  probability  of  coming  from  similar 
materials),  the  number  of  combinations  is  still  tremendous. 

A  practical  approach  is  to  treat  each  individual 
return  separately.  This  has  tremendous  practical  appeal 
because  of  reduced  calculation  requirements  and  has  less 
dependency  on  the  a  priori  probability  distributions. 

The  objective  of  processing  a  single  point,  for 
the  purpose  of  bistable  (bilevel) ,  grey  level  or  color 
displays,  or  even  for  coding  for  subsequent  retrieval, 
can  be  summarized  as  M-ary  detection.  In  other  words,  we 
want  to  assign  that  particular  point  to  one  of  the  M 
hypotheses  representing  a  specific  material.  Again  this 
technique  leads  to  the  well-known  Bayes  or  Neyraan-Person 
criterion.  Alternatively,  this  M-hypotheses  can  be 
expressed  in  terms  of  choosing  the  maximum  a  posteriori 
probability.  Without  loss  of  generality,  if  a  continuum 
of  grey  level  is  sought,  then  one  desires  an  MAP  (Maximum 
A  Posteriori)  estimate. 

As  mentioned  in  earlier  sections,  the  specular 
and  scatterer  parameters  vary  drastically.  Since  these 
underlying  parameters  are  unknowns  and  at  any  point 
different  probability  density  functions  can  occur,  it  is 
difficult  to  implement  the  aforementioned  techniques. 


Then,  it  is  desired  to  implement  a  robust  algorithm  with 
similar  signal  processing  approaches. 
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Without  loss  of  generality,  this  estimation  or  M-ary 
detection  problem  can  be  stated  as  assigning  the  received 
echo  s(*)  to  [o,lJ.  If  discrete  levels  are  desired,  such 
as  bilevel  displays,  computer  memory  or  color  display,  the 
resultant  assignment  is  further  processed  with  a  quantizer. 


where : 

s(*)  «R+  (positive  real  numbers), 
v ( • )  « [ 0 , l] ,  and 

V>  ‘fVo'Vl . M  ' 


I 


; 


64 


If  video  v(-)  is  used  for  grey  level  display,  then  without 
loss  of  generality,  let  0  be  dark  and  1  be  the  brightest 
level.  Let  VQ  be  dark  and  Vq  be  bright  also.  This  is 
not  the  only  way  to  assign  the  symbols,  especially  for 
color  displays.  It  is  beneficial  to  have  an  interactive 
system  where  an  operator  can  dynamically  change  the  color 
assignments  to  optimize  the  interpretation.  An  obvious 
assignment  for  color  display  is  to  let  red  represent  blood. 
This  grey  level>color  assignment  problem  is  not  essential 
to  the  dynamic  range  compression  algorithm  in  question. 

Now  we  can  delineate  the  approach  as  an  estimator 
which  has  the  form: 


s(-) 

+ 

Limiter 

v(-) 

i 

f 

f- 

E  o ,  i  ] 

g(') 

b( - ) 

L 

or 


vO)  =  tg(-)  s(-)  -  b  (  • )  J  L 

where  [  ]L  symbolizes  the  limiter  operation. 

Since  s(*)*R+,  the  magnitude  of  the  received  echo, 
the  effect  of  the  bias  level  b(-)  is  to  assign  s ( •  )<  [  0 -y] 
to  the  symbol  0.  The  effect  of  g(*)  is  to  assign 
s(  •  )<[  b(  • )  ,  b(  •  )  +  lJ  to  (0,1).  For  s(  •  )>fb  (•)  +  !?,  the  symbol 

g FT"  v  g(-)  ' 


1  is  assigned.  This  corresponds  to  the  dark  and 
saturation  levels  of  display  systems. 

To  repeat,  we  want  to  assign  the  point 
value  s(*)  to  a  specific  symbol  to  represent  a  certain 
type  of  biological  material.  This  simple  processing 
structure  will  accomplish  this  goal.  What  remains  is  to 
choose  the  bias  parameter  b(-)  and  the  gain  parameter 
g(’).  These  parameters,  of  course,  will  be  related  to 
the  underlying  probability  distributions. 

3.2.4  Bias  Parameter 

The  bias  parameter,  b(-)>  is  used  to  group  low 
level  echo  returns  in  one  common  symbol.  For  grey  level 
displays,  this  should  be  the  dark  level.  For  a  color 
display,  perhaps  red  is  appropriate.  In  echocardiography, 
these  low  level  returns  are  generally  from  a  volume  of 
blood  cells.  This  fact  might  not  be  true  for  other  medical 
ultrasonic  applications  such  as  obstetrical  abdominal  scans 
for  expectant  mothers,  where  the  lowest  returns  will 
probably  be  from  water.  But  nontheless  these  low  levels 
represent  returns  from  one  particular  material  and  are  of 
no  interest  for  medical  diagnosis.  In  radar,  these  returns 
are  generally  called  clutter,  where  one  possible  source 
might  be  returns  from  the  surface  of  the  ocean.  Then, 
this  bias  level  should  be  a  function  of  a  scatterer  para¬ 
meter  similar  to  that  from  blood  cells. 

As  mentioned  in  Section  2.3.1,  the  echo  from  a 
volume  of  blood  has  a  Rayleigh  distribution.  Furthermore, 
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we  saw  that  in  Section  2.3.2,  the  scattering  parameter 
varies  spatially  and  temporarily  due  to  attenuation,  wave 
divergence,  etc.  Since  these  parameters  are  unknowns  and 
we  seek  a  robust  algorithm,  this  bias  parameter  will  be 
derived  from  the  data  itself.  In  echocardiography,  failure 
to  detect  a  structure  (especially  valves  in  a  volume  of 
blood)  are  "costlier”  than  false  alarms.  This  is  especially 
true  for  an  image  because  the  probability  of  false  alarms 
over  a  small  neighborhood  is  highly  unlikely  and  thus  if 
a  false  alarm  occurs,  it  appears  at  an  isolated  point  and 
the  human  interpeter  can  easily  recognize  its  insignificance. 
This  couples  nicely  with  the  proposed  algorithm  since  this 
helps  in  eliminating  some  of  the  effects  of  acoustic  energy 
attenuation . 

The  following  is  an  examination  of  the  character¬ 
istics  of  the  bias  parameters  under  different  signal  para¬ 
meters.  Let 

PM  =  probability  of  missing  a  significant  structure 
and  M 

P  =  probability  of  false  alarm  (or  displaying  blood 
F 

volume)  . 

The  significant  structure  can  either  be  a  distinct  specular 
component  or  a  backscattering  from  tissue  layers.  Without 
loss  of  generality,  the  probability  density  function  for 
either  case  (H^)  is 
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Ps/h/s/H1> 


exp 


(s2*ni2) 

~2^J~ 


Sr|  -I 

V — &  ’ 

°1 


and  for  backscattering  from  the  blood  volume  (H^) , 


ps/Hrt^s/Ho)  =  2 

o  a 

o 


exp 


2o 


U(s) 


Since  the  bias  parameter,  b(-)=b,  operates  on  each  single 
point,  s,  the  probabilities  PM  and  Pp,  can  be  expressed  as 
follows : 


and 


•M 


(s/H1)ds 


(s/HQ) ds 


Solving  for  Pp, 


ds 
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One  can  fix  the  probability  of  false  alarm  P_ 

F 

by  choosing  the  threshold  b  proportional  to  o  .  In  this 

o 

algorithm, aQ  will  be  estimated.  Let 
b  =  Koq  , 


then 


Since  an  estimate  of  o  will  be  used  in  the 

o 

algorithm,  it  is  important  to  examine  the  sensitivity  of 


Figures  3.5  and  3.6  show  PF  vs .  k  and  vs.  k, 
respectively.  As  shown  on  the  graphs,  if  K>1,  one  can 
choose  K  to  give  relatively  good  performance.  But  this 
performance  is  in  terms  of  false  alarms  only,  it  is 
important  to  know  k  affects  missed  targets.  Since  the 
specular  or  scatterer  components  from  tissues  and  inter¬ 
faces  vary  widely,  P^  cannot  be  specified  explicitly.  We 
will  examine  how  P^  changes  as  the  parameter  changes. 
Recall  that  for  a  significant  target,  we  have: 
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(s2+n12) 


Ps/H,  <S/H1)  ‘  ^2  e  201  D(s> 

1  a,  a. 


Since  this  is  a  point-by-point  detection  problem  and  the 
threshold  value  is  b,  the  probability  of  a  missed  target, 


Ps/H.(s/Hl,dS 


00 

•  (  Ps/H  <S/H1> 

« J  ^  n 


=  1  -  /  —  e 

J  °1 


(s2+t^) 


20 1  I  (— j)  ds  • 


The  integral  cannot  be  evaluated  analytically. 
The  normalized  function  was  tabulated  by  Marcum  in  terms 
of  a  function  commonly  called  Marcum's  Q  function: 


C  2  2 

Q(a,B)  =  I  Z  exp(-  Z ■  -)  lo(0!z)dz  . 
B 

Therefore,  P  can  be  expressed  as 


PM  =  1  "  Q<-§»b)  • 
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If 


b 


Koo 


> 


then 


PM  -  1  -  «(~2  •  K°o>  . 

01 

which  can  be  expressed  as 

PM  =  1  ~  Q  <<~J2  ’  °o  /_2  -*n  PF  )  ' 
1 

Let 


Figure  3.7  shows  P„  for  various  normalized  values  of  y 

M 

and  | . 

Within  the  scanning  volume,  there  is  always  a 
region  of  blood.  Therefore,  one  can  estimate  oQ  at  that 
particular  region.  Recall  that  due  to  attenuation,  beam 
divergence  and  other  effects,  aQ  will  not  be  constant 
throughout  the  scanning  volume.  One  must  compensate  for 
this  effect  if  b=KoQ  is  to  be  used  throughout  the  whole 
volume.  Since  the  attenuation  effects  cannot  be  estimated, 
it  is  logical  to  compensate  for  the  beam  divergence.  This 
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Probability  of  Missed  Target 
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has  an  effect  of  using  a  lowered  threshold;  but  this  is 
desirable  because  of  reasons  mentioned  earlier. 

Estimation  of  oQ  depends  on  implementation 
techniques.  For  example,  real-time  processing  will 
force  an  analog  approach  rather  than  the  digital  approach 
used  for  nonreal-time  processing.  Also,  the  optimization 
of  hardware  will  dictate  the  estimation  technique.  Since 
echo  signals  from  this  volume  of  blood  has  a  Rayleigh 
distribution,  the  estimate  can  be  the  sample  mean  or  from 
second  order  statistics,  which  follows  from, 

-/¥  °o 


and 

Ep/Koi  -  2oo2  . 

For  calculating  g(*),  we  need  the  higher  order  statistics, 

therefore,  from  the  hardware  standpoint,  it  is  better  to 

estimate  oQ  with  a  mean  square  value.  Since  the  blood 

scattering  occur  over  a  volume  (spatial  and/or  temporal) 

the  estimate  can  best  be  formed  over  the  whole  volume 

depending  on  the  dimensionality  of  the  scan.  The  large 

number  of  samples  needed  to  make  the  estimates  tends  to 

decrease  the  error  (proportional  to  — ) .  For  computer 

y/N 

approach,  a  running  estimate  over  the  entire  volume  and  a 
comparison  to  stored  values  of  the  minimum  will  be  used. 
The  "box  car"  estimate,  o  , 


can  be  expressed  as 
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A. 
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Nt  N  n  nt 
i+~T  k+_5~  1+T“ 


i _ . _ y  y  y  y  ■ 

2niWl  4  ^nt  Z"ni  n. 

w-i-2^  x-j ^  y=-k— j-  z=l— ^ 


s  (w,x,y ,z) 


1/2 


where  Nj ,  Nj,  NR,  denote  the  number  of  samples  over 
the  respective  dimensions,  i,  j,  k,  1.  Obviously  the 
total  number  of  samples  is  equal  to  NTxNTxN„.xNT  .  For 
practical  purposes,  and  in  the  case  of  snapshots,  the 
i-axis,  or  time,  is  not  used,  and  NT=NT=NV  since  no 
obvious  advantage  can  result  if  one  particular  axis  is 
favored.  A  trade-off  in  using  a  larger  volume  is  between 
computation  time  and  the  geometric  dimensions  of  the  heart 
chamber.  A  good  reliable  estimate  can  be  calculated.  For 
real  time  application  applications,  averaging  over  range 
can  be  accomplished  by  utilizing  an  analog  filter;  this  will 
enable  a  lower  sampling  rate  for  this  channel  for  the 
purpose  of  computing  b  and  g.  Figure  3.8  illustrates  such 
a  scheme.  Here,  the  computer  or  digital  filter  is  used  to 

A 

store  adjacent  lines  and  to  compute  oq  over  the  other 

dimensions.  Now,  the  minimum  value  of  a  can  be  found  for 

o 

an  estimate  of  for  blood.  This  value  is  valid  at  that 
o 

particular  volume  only,  and  as  shown  in  Section  2.3.2,  we 
see  that  oq  is  a  spatial  and  temporal  function  due  to  the 
beam  divergence,  attenuation,  etc.  An  example  of  the  bias 
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value  adjusted  for  the  beam  divergence  is  illustrated  in 
Figure  3.9. 

3.2.5  Gain  Parameter 

The  effect  of  the  bias  parameter,  as  described 
previously,  is  to  assign  all  low  level  returns  from  blood 
to  one  symbol.  On  the  other  hand,  the  gain  parameter  has 
two  major  effects.  Analogous  to  the  bias  parameter,  the 
gain  parameter  tends  to  assign  all  "high"  returns  from 
specular  reflectors  to  a  common  symbol.  Also,  depending 
on  the  instantaneous  gain  value,  the  gain  parameter  can  be 
used  to  assign  a  specific  code.  In  other  words,  besides 
the  saturation  effects,  the  gain  parameter  couples  with  the 
quantization  scheme. 

The  major  objective  of  the  gain  parameter  is  to 
maximize  the  information  content  of  the  dynamic  range- 
compressed  image.  This  is  similar  to  a  source  coding 
problem.  Due  to  the  subjective  nature  of  a  fidelity 
criterion  for  images,  a  classical  source  coding  criterion 
approach  is  not  taken  here.  Another  reason  for  not  using 
the  source  coding  approach  is  the  highly  nonstationary  nature 
of  the  echo  returns.  First,  we  will  decide  which  approach 
should  be  taken  if  a  priori  statistical  information  is 
available  for  the  volume  scanned. 

As  described  in  the  previous  section,  the  statistical 
model  of  the  ultrasonically  scanned  volume  consists  of 
Rayleigh  distributions  due  to  scatterers  and  Rician  dis¬ 
tributions  due  to  specular  reflectors.  Then,  setting  the 


Figure  3.9  Example  of  Beam  Divergence  Compensation 


79 


gain  parameter  involves  normalizing  the  distribution  so 
that  a  satisfactory  maximum  value  will  be  equal  to  the 
saturation  level  of  the  encoder,  for  example,  an  analog- 
to-digital  converter.  Then,  if  we  consider  that  a 
percent  of  high  level  echo  returns  can  be  grouped  under 
one  symbol,  the  saturation  level,  one  easily  can  evaluate 
the  gain  parameter  at  each  point.  We  consider  only  a 
Rician  distribution  only  for  the  following  discussion 
(Rayleigh  distribution  is  just  a  special  case) .  The 
probability  of  the  echo  signal  being  higher  than  a  value 


where  Q  is  the  Marcum's  Q  function.  Clearly, 


M 


=  f(a, 


V 


2 

where  r\  and  a  are  from  the  normalized  distribution 
'o  o 

function.  Now, for  an  arbitrary  point  in  the  scanning 

2 

volume  with  r\  and  o  ,  the  saturation  level,  sM,  is: 


sm  =  f(a’  tl»  °2)  * 


The  gain  parameter  can  be  found  easily  as 
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8  '  Vsm  ' 

Using  this  gain  parameter,  the  echo  returns  from  significant 
reflectors,  either  specular  or  scatter  components,  tend  to 
form  a  stationary  process.  The  exception  occurs  at  regions 
of  blood  scatterers.  With  the  application  of  the  bias 
parameter,  the  blood  regions  will  have  zero  intensity. 
Therefore,  the  total  effect  is  that  the  significant  cardiac 
structures  are  evenly  illuminated  by  the  acoustic  wave. 

Besides  the  unavailability  of  the  a  priori  statistics 
other  problems  are  present  in  implementing  this  gain  parameter. 
For  abrupt  changes  such  as  from  endocardium  to  blood  volume, 
the  gain  tends  to  switch  from  a  low  value  to  a  high  setting. 

In  practical  hardware  implementation,  this  abrupt  change 
will  tend  to  generate  an  artifact,  a  false  return.  A  more 
gradual  change  in  gain  can  combat  this  problem.  Low  pass 
filtering  and  using  continuous  gain  profiles  are  ways  to 
implement  this  gradual  change  in  gain  profile.  Besides  the 
need  for  a  priori  statistics,  to  implement  this  scheme,  one 
must  also  know  where  the  blood  volume  occurs.  This  does 
not  offer  a  practical  implementation  method  which  can 
be  applied  to  the  whole  scanning  volume.  But  all  is  not 
lost  since  the  information  from  the  bias  parameter  offers 
a  clue  to  the  location  of  blood  volume;  then,  the  gain  para¬ 
meter  can  be  set  with  the  bias  parameter  as  one  of  the 
dependent  variables.  These  points  are  taken  into  account 
in  the  algorithm  proposed. 
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The  gain  parameter  can  be  calculated  with  the 
following  steps: 

* 

1.  Use  the  "box  car"  o  estimate  from  the 

n 

bias  parameter  calculation.  This  saves 
computational  requirements. 

2.  Based  on  the  ratio  a  /a  .  ,  where 

o  o ,  mm  . 

A  A 

a  denotes  min.(a  )  over  the 

o , min  o 

scanning  volume,  form  a  multiplier  k  from 
a  nonlinear  function,  fNL- 

3.  Form  the  multiplier  g  =  s  /ko  ,  where 

IQHX  O 

s  is  the  saturation  level, 
max 

The  nonlinear  function,  fv  is  such  that 

r'lL 


A  ^ 

f  XTT  (o  /o  .  =  1)  =  k 

NL  o  o  y min  1 


and 


fNL<S</Vmin-?  “>  *  k«  ' 

with  k  >k  .  For  example,  a  piecewise  linear  function 

O  CO 

illustrated  in  Figure  3.10  is  effective. 

k  is  chosen  such  that  the  saturation  effect  is 
x 

satisfactory.  An  estimate  of  what  k^  should  be  will  be 

covered  later  in  this  subsection.  If  kQ  is  chosen  greater 

than  k  ,  the  effect  is  such  that  the  resultant  brightness 
00 

level  of  the  display  is  reduced  for  certain  regions  where 
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'o  W'O  The  regions  where  this  will  occur  are  within 

U  W  y  IU  J.  XI 

the  blood  volume  itself  since  o<&o  .  .  in  that  case, 

o  o,mm  ’ 

there  are  two  reasons  why  k  <k  .  First,  S  will  be  on  the 
low  side  at  the  isolated  structure  due  to  "box  car" 
averaging.  The  second  reason  is  that  this  isolated 
structure  occurs  in  a  volume  of  darkness,  that  is,  blood. 

In  order  to  reduce  artifacts,  the  brightness  level  should 
be  subdued.  Of  course,  kQ  should  be  operator  interactive. 
This  does  not  change  the  adaptive  nature  of  the  algorithm; 
it  merely  gives  the  operator  a  choice  for  visualizing 
isolated  structures  better  at  the  possible  expense  of 
increased  artifacts. 

Another  way  to  interpret  this  algorithm  is  to 
consider  oq  as  a  root  mean  square  power  estimate. 
Normalization  using  this  power  estimate  is  a  common 
practice.  The  nonlinearity  is  introduced  to  reduce 
brightness  somewhat  in  the  blood  region. 

The  question  remains  on  how  this  algorithm 
works.  A  statistical  analysis  of  how  this  receiver  behaves 
is  presented  here.  A  simplification  is  used  to  illustrate 
the  receiver  behavior  that  is,  only  a  scattering  component 
is  considered.  The  analysis  applies  since,  over  a  local 
volume,  the  actual  distribution  function  is  a  conglomerate 
of  different  Rayleigh  and  Rician  distributions,  and  thus, 
the  "sharper"  distribution  is  "broadened",  as  discussed 
previously . 

Let  us  consider  a  small  volume  where  the  distribution 
is  Rayleigh,  or  the  returns  are  mainly  from  a  heterogeneous 
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and/or  layered  region.  Also,  oQ  for  this  region  is  much 
greater  than  o  .  for  the  scanning  volume.  Then,  the 

U  f  111  X  11 

gain  parameter : 

s 

„  max 

s  =  T — “  • 

k  a 
oo  o 

Then,  the  echo  returns  in  this  region  that  exceed  the 
saturation  level  are  such  that 


gs  >  s 


max 


or 


or 


max 

.  A 

k  a 
oo  o 


>  s 


max 


A. 

s  >  k  o 
oo 


o 


Since 


9 


the  probability  or  portion  of  the  echo  in  the  saturation 
level  is: 
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Suppose  aQ  =  oQ,  and  solving  for  : 

koc  “  S~2J  n  PSAT  • 

For  example,  if  5%  saturation  is  satisfactory,  then 

k^  =  /-2  JLn(  .051)  =  2.45.  Figure  3.11  shows  a  plot  of 

P.,_  versus  k  .  It  is  important  to  examine  the  effects  of 
SAT  oo 

o  estimation  error.  Let 
o 

°o  "  ao(1+7)- 

where  7  =  fractional  error  is  estimating  aQ. 


Then,  the  actual  probability  of  saturation, 


1)  Ap  deviates  greatly  as  y  increases  since 

k  >1, 

00  ’ 

2)  Underestimating  aQ  implies  a  greater 
portion  of  the  signal  in  the  saturation 


level  and  vice  versa. 
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Since  aQ  is  estimated  from  box  car  averaging,  there  are 
basically  two  types  of  regions  which  will  give  /aQ 
estimation  errors.  They  are: 

1.  Edge  effects.  When  the  averaging  volume 
is  taken  over  a  region  where  the  endo¬ 
cardium  meets  the  blood  chamber,  an  over¬ 
estimation  will  occur  on  the  blood  side 
and  an  underestimation  will  occur  on  the 
cardiac  wall  side.  The  overestimate 

will  cause  the  initial  blood  region  to  take 
on  a  darker  value,  whereas  the  under¬ 
estimate  will  cause  the  edge  of  the  cardiac 
wall  to  take  on  a  brighter  value.  This  is 
essentially  an  edge  enhancement  effect. 

This  is  equivalent  to  an  additive  differ¬ 
ential  component.  Fortunately,  this  occurs 
as  a  positive  effect. 

2.  Isolated  points.  This  can  either  be  an 
isolated  reflector  in  a  blood  volume  such 
as  a  valve  or  a  small  fluid  volume  such  as 
the  epicardian.  In  the  former  case,  we  see 
that  is  definitely  an  underestimation 
since  the  box  care  region  is  larger  than  the 
isolated  point.  As  a  matter  of  fact,  this 
estimation  error  will  be  worse  than  the  edge 
effect.  With  the  nonlinearity  introduced, 
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this  effect  can  be  reduced  somewhat.  But  the 
overall  effect  is  that,  for  the  isolated  point, 
the  brightness  level  will  be  high.  This  is 
highly  desirable.  For  the  case  of  the  small 
fluid  region,  the  effect  is  exactly  opposite: 
we  have  a  darker  region.  Again,  this  is  a 
desirable  result. 

Then,  from  this  analysis,  an  extra  benefit  of  the 
dynamic  range  compression  algorithm  is  image  enhancement. 

3.2.6  Quantization 

This  discussion  pertains  to  echocardiographic 

imaging  systems  where  the  echo  returns  are  sampled  and 

stored  by  digital  means  for  subsequent  image  processing 

and  viewing.  This  can  further  be  separated  into  systems 

where  image  processing  techniques  are  used  before  coding 

into  brightness  levels  or  pseudo-color  assignments  and  to 

systems  where  the  digital  memory  performs  as  an  image  storage 

medium  or  scan  converter  only  with  no  numerical  processing 

of  the  echo  data.  The  difference  between  the  two  cases  is 

that,  in  the  former,  the  quantizer  has  a  narrower  dynamic 

range  than  in  the  latter.  For  example,  if  the  dynamic 

range  compressed  echo  data,  as  described  in  this  chapter, 

is  sampled  at  2  MHz  and  quantized  to  eight  bits,  then  the 

dynamic  range  of  the  compressed  signal  must  be  less  than  48 

dB.  This  can  easily  be  accomplished  since  for  small  range 

(69) 

segments,  the  dynamic  range  is  of  the  order  of  30-50  dB.  ' 
The  dynamic  range  compression  algorithm,  acting  as  a  local 
normalization  scheme,  can  easily  compress  the  50-70  dB 
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variation  due  to  attenuation.  For  this  case,  where  the 
data  is  further  processed,  the  quantization  scheme  is  well 
described  in  the  literature  mentioned  in  Section  3.1. 

For  practical  grey  scale  displays,  where  the  dynamic  range 
is  limited  to  15  to  20  dB's,  or  for  pseudo-color  displays 
where  the  colors  available  are  limited,  the  immediate 
quantization  of  the  compressed  data  should  be  "optimized.” 

Of  course,  the  actual  processing  can  be  separated  into  two 
stages.  First,  the  dynamic  range  compressed  data  is 
quantized  into  eight  bits.  Then  the  "quantization"  scheme 
can  be  implemented  again  with  the  computer.  But  if  the 
eight-bit  data  are  available  in  the  computer  already,  other 
image  "enhancement"  schemes  such  as  histogram  equalization 
can  be  applied.  But  the  disadvantage  is  the  processing 
time  required.  For  real-time  systems,  this  is  not  feasible. 
Therefore,  this  section  addresses  to  the  problem  of 
quantizing  directly  to  display  level  dynamic  ranges.  Bear 
in  mind  that  the  resultant  code,  the  quantized  value,  is 
not  the  video  voltage,  but  rather  a  grey  level.  In  other 
words,  a  gamma  amplifier,  or  an  antilog  amplifier,  is 
incorporated  into  the  grey  level  display  to  compensate 
human  visual  responses. 

In  order  to  optimize  the  quantization  scheme,  one 
must  first  choose  a  reasonable  criterion.  If  the  purpose 
is  data  processing,  a  mean  square  cost  is  a  reasonable 
choice.  Another  common  choice  is  that  constant  percentage 
error  is  desired.  This  leads  to  a  logarithmic  quantization 
scheme  commonly  used  for  image  processing.  For  echocardio- 
graphic  imaging,  this  scheme  might  not  be  statistically 
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optimal.  As  a  matter  of  fact,  such  images  usually 
improve  with  local  histogram  equalization.  Then,  a 
reasonable  choice  of  criterion  is  the  average  information 
contained  in  the  quantized  data.  Then  for  the  quantizer, 


s 


Quantizer 


let  input  s,  s  €  (°»s^x)  and  output  q,  q  =  i  if 

d  <  s  <  d  for  i=0 , 1 , . . . , N-l .  Usually  N=2B;  for 
i-1  1 

example,  a  16-level  grey  scale  display,  N=16  and  B=4. 

The  set  of  values  are  decision  levels  which  have  to 

be  decided. 

With  the  information  theoretic  approach,  then 

we  want  to  maximize  the  mutual  average  information  between 

the  input  and  the  output.  This  is  an  application  of  a 

discrete  memoryless  channel  (DMC)  with  the  source  alphabets 

/ 13) 

extended  to  an  infinite  set.  It  is  well  known  '  that  for 
a  symmetric  discrete  memoryless  channel,  full-channel, 
full  capacity  can  be  obtained  by  using  the  inputs  with 
equal  probability.  By  partitioning  the  set  of  infinite 
source  alphabets  into  subsets  equal  to  the  number  of 
quantization  levels,  it  follows  that  the  full-channel 
capacity  is  used  or  the  average  mutual  information  is 
maximized.  Therefore,  the  decision  level  should  be  chosen 


such  that  the  probability  of  the  quantization  input  falling 
into  that  interval  is  equal  to  1/N,  the  inverse  of  number 
of  quantization  levels. 
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A  direct  derivation  of  such  a  result  are  as  follows: 


Let 


p_(s)  =  input  probability  density  function, 


Pg(qn)  =*  P(q-n)  =*  output  probability,  and 
I(s,qn)  =  average  mutual  information. 


Then , 


I(s,qn) 


N  r° 

■S  / 1 

n=l  J 
-oo 


PQ(qn)pS  Q(S  qn)log 


PS  Q(S  V 
p(s) 


ds . 


Note  that,  if  s  is  quantized  according  to  decision  levels 
dn  ,  then ; 

d 

n 

W  -  f  ps<t)dt 

d  , 
n- 1 

and 


P(s/qn) 


Pg(s)  l  U(s-dn_1)-U(s-dn)J 


Ps(t)dt 


Now  , 


PTTvT  PS(S)  ^D(s-dn-l)-U(3-dn)f 


n  ?  r 

'll  I  PQ<V  p-TTTPS(s)[D<s'dn-l)-U(s-dn> 

n=l  J  W  n  * 

-x 


*  log 


WJ  ps(s)  [«<s-dn.1)-U(s-dn)] 


pS(s; 


N  y.n 

I(s,qn)  =  “2^  J  Ps(s)log  PQ(qn)ds 


d 

N  -n 

S ./  ■ 


J  Pg (s) log  /  Pg(t)dt  ds  . 
dn-l  dn-i 


Maximizing  I(s,qn>  over  > 


n 

St-  I(s>qn)  “  -PS(dn)log  f  Ps(t)dt+PS(di> 

0  A  J 


+  p„(d_)  log  /  p_ ( t) dt-p0 (d . )  , 
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Setting 


Or 


aT-  KS’V  -  0 

n  dn+l 


log 


■/P 


p  (t)dt 
d  7  s 


J 


n 


Ps(t)dt 


/ 


n-1 


n-1 


Pg(t)dt  = 


n+1 

/  Ps' 


which  is  the  desired  result.  Since 


i  =  1,2, . . . ,N 


Pg(t)dt  , 


N  N 

T"!  J  Pg(t)dt  =  y  Pg(t)dt  =  1  , 

dn-l 


then  , 


d 

/ 


ps(t)dt  =  £  • 

n-1 

One  way  to  implement  this  is  to  use  an  instantaneous  trans¬ 
formation,  T(s)  ; 

s 

w  =  T(s)  =  j  Pg(t)dt  , 

-oo 

and  feed  w  into  a  uniform  quantizer.  Note  that  w«Lo,l]  since 
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it  is  a  cumulative  distribution  function,  and  for 
wc  [  0 , 1) ]  : 


pw(w) 


I®/ 


_1(oo)) 

X“1 

Ps(t) dt 


“00 


Pg  (S=T_1  (o>)  ) 
Pg (S=T_1 (w) ) 


Then, the  transformation  maps  the  input  into  a  uniformly 
distributed  output  over  Co,l].  Since  the  uniform 
quantizer  has  £dn^  =  ^0,^,^,...,!?  ,  then; 


n 

/ 


Pg ( t ) dt 


n-1 


Pw(t)dt  = 


1 

N 


It  is  interesting  that  this  result  is  the  same  as  quanti- 

/  j  o  \ 

zation  with  a  minimum  mean  square  error  distortion  criterion. 

The  results  were  derived  in  this  case  avoiding  the  definition 
of  an  output  symbol  and  the  average  mean  square  error 
criterion.  For  the  transformation  of  raw  data  into  image 
data,  this  has  a  more  intuitive  appeal.  As  a  consequence  of 
this  quantization  scheme,  if  the  data  are  used  for  numerical 
processing,  the  quantization  scheme  is  still  optimal  in  the 
least  mean  square  sense. 

For  application  of  this  maximum  entropy  quantization 
scheme  to  this  dynamic  range  compression  algorithm,  a  non¬ 
linear  amplifier  should  be  used  in  front  of  the  analog-to- 
digital  converter  (a  uniform  quantizer) .  Since  the  gain 


96 


parameter  "tends"  to  normalize  the  data,  a  time-invariant 
nonlinear  amplifier  will  suffice  due  to  the  pseudo¬ 
stationary  nature  of  the  compressed  signal.  The  nonlinear 
function  can  be  derived  by  assuming  the  input  distribution 
to  be  Rayleigh  (over  a  small  neighborhood) .  One  reason  is 
that  there  is  a  lack  of  information  for  the  Rician  distri¬ 
bution  due  to  specular  reflectors.  With  the  normalization 
scheme  applied,  the  specular  return  tends  to  be  at  the  high 

i 

levels  (bright  levels  for  a  grey  scale  display);  thus,  the 
quality  of  the  image  itself  is  not  sacrificed.  Then,  for  a 
Rayleigh  distribution: 


Pg(s)  =  — S~2  expL  -  — — g-J  U(s) 
°o  °o 


Then  the  transformation  desired  is: 


T(s)  = 


C 


— 2  exP^  ”  dt 

a  o 

o  °  2  ° 

s 


=  1  -  e 


% 


The  input-output  characteristic  is  shown  in  Figure  3.12. 

Now  since  the  input  to  a  "real  world”  amplifier  has 

a  limit,  which  was  in  the  last  subsection  called  S_„„,  this 

max 

should  be  set  slightly  above  the  d^_^  decision  level. 

Since  the  output  is  uniformly  distributed,  and  the  quantizer 
has  uniform  decision  levels: 


d 


N-l 


N-l 

N 


Then 


S 


max 


Vi +  A 


9 


where  A  -  small  quantity  choose  by  the  designer 
and  has  no  significance,  and: 


.  _  N-l 

*N-1  =  N 


N-l 
Z 


1  - 


2  a 

e  o 


Solving  for  tN-1, 


Vl  -  /2°o  N  ' 


and  therefore, 


Smax  N  +  A 


Or  if  S  is  fixed  by  the  equipment,  then: 
max 


n  2  _  1  c;2 

Joo  ■  Jn^  “max  * 


Then, the  transformation  required  is  ; 


where  b  =  bias,  and  K  =  normalization  constant.  This 
is  due  to  the  fact  that  the  signal  has  been  biased 
to  remove  the  blood  component.  Since  b  is  usually 
much  smaller  than  aQ,  the  assumption  that  b=0  is  of  no 
great  consequence  and  is  much  easier  to  implement. 

Another  note  of  importance  is  that  one  symbol, 
say  symbol  0,  has  been  used  for  blood  level  already. 

That  is,  if  s=0,  then  this  represents  blood.  If  s/0 , 
that  means  the  signal  is  of  interest,  and  this  should  be 
represented  by  a  symbol  other  than  0.  For  example,  if  a 
quantizer  has  N  quantization  levels  on  the  interval  [ 0 , lJ , 
we  should  consider  quantizing  s  into  N-l  symbols  other 
than  the  symbol  0.  This  can  be  accomplished  by  adding 
an  offset  and  by  multiplying  by  a  scale  factor.  Then, 
this  second  transformation  is: 


T2(s)  “  ^offset  +  Kmult 


s, 


100 


where 


and 


Also, 


Kof fset  N 

TT  N-l 

mult  N  ' 

the  first  transformation  should  be  modified  to: 
_  In(N-l)  g2 
Tx(s)  -  1  -  e  S  max 


be : 


Combining  T^,  Tg ,  the  overall  transformation  should 

inN-1  _2 

'  r2 —  s 

T(s)  -  1  -  e 


Then,  if  T(s)  is  fed  into  a  uniform  quantizer  with  decision 
levels  (0,  i,  ^,  .  .  .  ,  ,  and  with  corresponding  symbols 

(0,  1,  2,...,  N-l),  the  probability  distribution  of  symbols 
around  a  local  neighborhood  will  tend  to  be  uniform,  which 
is  the  intention  of  the  dynamic  range  compression  algorithm. 

3.2.7  Implementation  Techniques 

The  dynamic  range  compression  algorithm  is  basically 
a  local  normalization  and  information  content  optimization 
technique  with  inherent  real-time  application  capabilities. 
The  actual  implementation  of  such  a  technique  to  an  echo- 
cardiographic  imaging  system  depends  on  the  system  itself. 
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It  is  a  feedforward  processor  with  its  advantages  of 
possible  infinite  system  gain  and  avoidance  of  stability 
problems.  Although  conceptually  it  is  a  two-channel 
device,  it  can  be  implemented  as  a  single-channel  device 
with  a  gain  parameter  acquisition  scan  to  derive  the 
control  signals  rather  than  using  a  separate  channel.  In 
this  subsection,  some  real  time  processor  architectures  are 
suggested  for  different  scanning  systems.  Some  basic 
timing  considerations  are  examined  to  show  their  feasibility. 

Figure  3.13  shows  a  possible  configuration  for  real¬ 
time  processing.  The  receiver  has  two  channels:  one  channel 
for  video  data  and  one  channel  for  deriving  gain  and  bias 
parameters.  The  dynamic  range-compressed  video  data  can  go 
directly  to  the  display  unit  or  the  quantizer  for  computer 
storage  or  processing.  A  separate  detector  can  be  used  in 
the  processor  channel  so  that  fidelity  of  the  video  data 
can  be  improved  by  detection  after  dynamic  range  compression 
on  the  video  channel.  The  log-amp  and  low-pass  filter 
combination  is  used  to  derive  the  power  statistic  within 
scan  line  segments.  The  low-amp  compresses  the  dynamic 
range  to  accommodate  the  dynamic  range  to  the  analog-to- 
digital  converter.  Here,  optimal  compression  and  quantization 
is  not  necessary  since  only  an  estimate  of  the  power 
statistics  is  required.  The  low-pass  filter  provides 
averaging  over  range  (or  equivalently  time)  and  thus,  the 
sampling  rate  requirement  for  the  A/D  converter  is  reduced. 

For  example,  if  sixteen  segments  of  gain  parameters  are  used, 


Figure  3.13  Two-Channel  Real-Time  Dynamic  Range  Compression 
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which  clinically  is  shown  to  be  adequate,  the  sixteen 
samples  required  over  the  200  nSec.  maximum  range  of  15 
cm  means  that  the  sampling  rate  is  80  KHz.  Therefore, 
a  separate  channel  without  the  higher  cost  MHz  sampling 
rate  A/D  converter  is  feasible  with  just  a  low  cost  A/D 
converter.  Since  the  conversion  time  for  the  A/D  converter 
is  much  longer,  it  is  possible  to  use  more  bits  to 
improve  dynamic  range.  After  the  power  statistic  data  has 
been  digitized,  it  is  stored  in  an  appropriate  memory  or 
registers  to  provide  averaging  in  the  other  dimensions. 

The  processor  can  be  dedicated  hardware  for  speed  or  a 
microprocessor  for  economy  depending  on  the  speed  require¬ 
ments.  The  processing  can  be  accomplished  between  segments 
samples  or  lines.  For  example,  M  mode  echocardiographic 
units  usually  have  a  PRF  (pulse  repetition  frequency)  of 
about  1  KHz,  since  the  data  interval  per  line  is  only  200 
nSec.,  about  800  (iSec.  are  available  before  the  next  pulse 
occurs.  An  economical  microprocessor  can  use  this  time  to 
derive  the  gain  and  bias  parameters  and  to  load  the 
registers  for  waveform  generation  with  no  severe  constraint 
on  software  considerations.  But  for  multidimensional 
sector  scanners,  the  timing  requirement  is  much  more  severe. 
Due  to  the  high  data  density  required,  PRF  is  usually  set 
at  a  maximum,  which  means  no  available  time  between  receiving 
time  and  the  next  transmit  pulse.  Therefore,  for  the  same 
sixteen-segment  example,  processing  must  be  accomplished 
within  the  time  for  segment  samples,  which  is  12.5  jiSec . 

For  a  simple  digital  filter  with  dedicated  hardware,  this 
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can  be  synthesized.  In  summary,  the  processor  calculates 
the  gain  and  bias  parameter  with  power  statistics  from  the 
ith  and  previous  frames  and  this  set  of  parameters  is 
applied  to  the  i+lth  frame. 

A  similar  application  of  this  dynamic  range  compression 
algorithm  can  be  applied  with  one  channel  only.  Figure  3.14 
shows  such  a  configuration.  This  is  especially  applicable 
to  a  snapshot  camera  system  since  only  one  frame  of  data 
is  required.  The  procedure  can  be  divided  into  two  scans. 

The  first  scan  acquires  the  power  statistics  for  calculating 
the  gain  and  bias  parameters.  The  second  scan  does  the 
actual  scan.  Since  the  physicians  are  interested  in  one 
particular  time  instant  of  the  cardiac  phase,  an  electro¬ 
cardiogram  is  used  to  synchronize  both  scans.  During  the 
first  scan,  zero  bias  and  a  norminal  multiplier  waveform  is 
used  to  compensate  for  the  gross  attenuation  effects.  The 
optional  logarithm  amplifier  can  then  be  switched  in  to 
ensure  the  dynamic  range  falls  within  the  range  of  the 
analog-to-digital  converter.  Again,  since  only  an  estimate 
of  the  power  profile  is  required,  the  logarithmic  compression 
is  sufficient.  The  applied  gain  profile  and  the  compressed 
echo  data  can  then  be  used  to  calculate  the  gain  and  bias 
parameters.  This  new  set  of  parameters  is  used  in  the 
second  scan  along  with  the  "optimal"  nonlinear  amplifier. 

One  advantage  of  this  configuration  is  that  the  method  and 
constants  used  to  calculate  the  gain  and  bias  parameters  can 
be  changed  easily  with  the  software  package  of  the  computer. 


Figure  3.14  One-Channel  Dynamic  Range  Compression 


Furthermore,  with  the  use  of  zero  bias  and  with  the  gain 
parameter  applied,  the  full  dynamic  range  data  can  be 
reconstructed.  This  is  especially  useful  for  research 
purposes  such  as  beam  width  deconvolution.  This  recon¬ 
structed  full  dynamic  range  data  definitely  offers  processing 
advantages  over  the  currently  used  data  with  arbitrary  and 
unknown  bias  and  time  varying  gains. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 

In  Chapter  II,  an  echocardiographic  signal  model 
is  presented.  In  this  chapter,  specific  examples  of  the 
echocardiographic  signal  model  are  described  together  with 
computer  simulation  results.  The  computer  results 
simulate  power  profiles  of  M-mode  echocardiographic  signals 
(from  stationary  and  dynamic  targets),  and  also,  three- 
dimensional  sector-scanned  snapshot  echocardiographic 
signals.  These  results  were  needed  to  check  dynamic 
range  compression  algorithms  in  two  respects.  First,  the 
signals  should  correspond  to  fairly  complex  geometries  to 
be  a  representative  test  of  the  algorithm.  (Because  of 
the  complexity  of  the  geometry  and  the  need  for  dynamics, 
it  was  impractical  to  set  up  in-vitro  experiments  to  gather 
repeatable  echocardiographic  signals.)  Second,  it  was 
not  economically  feasible  to  record  the  wide  dynamic 
range,  high  speed  raw  echocardiographic  signals  to  form 
a  sufficiently  comprehensive  data  base.  To  overcome 
these  obstacles,  these  simulated  echocardiographic  signals 
are  used  in  the  chapter  to  examine  the  adaptive  dynamic 
range  compression  algorithms. 

In  Chapter  III,  a  general  dynamic  range  compression 
applicable  to  multidimensional  echocardiographic  imaging  of 
any  dimension  is  presented.  Implementation  of  the  algorithm 
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for  one-dimensional  real  echo  data  is  shown  in  Appendix  A. 
In  this  chapter,  specific  implementation  methods  of  the 
algorithm  are  presented.  Due  to  the  nature  of  hardware 
and  of  the  economics  involved,  implementation  of  the 
adaptive  dynamic  range  compression  algorithm  takes  on 
different  forms  for  different  scanning  and  display 
techniques.  This  is,  of  course,  at  some  expense  to  the 
"performance"  of  the  algorithm.  At  this  stage  where  no 
viable  automatic  gain  control  is  available,  especially  for 
sector  scanners,  a  small  sacrifice  in  "performance"  is  not 
detrimental  to  the  objectives.  Besides,  no  objective  per¬ 
formance  index  for  echocardiographic  images  is  available. 
One  who  operates  an  imaging  system  must  be  aware  of  the 
system  and  its  behavior,  and  must  realize  that  the  image 
generated  is  a  signature  rather  than  a  facsimile  of  the 
cardiac  geometry.  Then,  if  the  dynamic  range  compression 
algorithm  maximizes  the  information  conveyed  to  the 
operator  without  serious  untractable  distortions,  the 
algorithm  accomplishes  its  design  objectives. 

In  this  chapter,  an  economical  M-mode  dynamic  range 
compression  technique  and  a  dynamic  range  compression 
technique  for  a  real-time  three-dimensional  sector-scanned 
snapshot  camera  are  presented.  Economy  in  the  M-mode 
implementation  is  accomplished  by  hybrid  processing.  The 
analog  integration,  and  digital  sampling  technique 
eliminates  the  need  for  high-cost,  high-speed  analog-to- 
digital  converters.  Thus,  this  technique  is  feasible  for 
augmenting  existing  M-mode  echocardiographic  imaging 
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systems.  For  the  three-dimensional  snapshot  camera  under 
exploratory  development  at  the  Applied  Research  Laboratory 
at  The  Pennsylvania  State  University,  digital  storage  is 
utilized.  Therefore,  the  main  problem  of  concern  is  the 
processing  time  required.  To  minimize  processing  time,  a 
table  look-up  technique  is  proposed.  An  initial  scan 
acquires  the  parameters  necessary  and  a  variable  look-up 
table  is  computed.  During  the  diagnostic  scan,  the  table 
can  then  be  dynamically  interrogated  and  a  resultant 
equivalent  dynamic  range  compression  applied.  Some  dis¬ 
cussions  of  how  these  techniques  behave  are  included. 

4. 1  Simulation  of  Echocardiography  Signal 

To  represent  typical  echocardiographic  signals 
encountered  in  clinical  applications,  the  echo  signals 
should  include  the  following:  1)  a  region  of  blood,  2) 
isolated  specular  component  within  the  blood  volume  to 
simulate  heart  valves,  and  3)  varied  angles  of  incidence. 

These  attributes  are  necessary  to  check-out 
saturation,  wash  out  and  transient  effects  of  dynamic  range 
compression  algorithms.  The  last  two  items  represent 
typical  problems  the  operator  faces  in  order  to  adjust  for 
a  usable  echocardiographic  image. 

An  algorithm  is  presented  here  to  generate  a  set  of 
echocardiographic  signals  given  a  specific  scanning  method 
are  the  specific  cardiac  structure.  M-mode  data  from 
stationary  targets  and  moving  targets  are  simulated  with  the 
algorithm.  Three-dimensional  sector-scanned  echocardiographic 
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data  are  also  presented.  These  data  are  used  in  the  evalu¬ 
ation  of  the  dynamic  range  compression  algorithm. 

4.1.1  Simulation  Algorithm 

Two  algorithms  are  used  to  simulate  echocardiographic 
signals.  First,  given  a  set  of  geometric  descriptions  of  the 
cardiac  structure,  with  the  associated  biological  parameters, 
the  specular  and  scatterer  components  and  attenuation  are 
obtained  as  a  function  of  range  for  each  transmit/receive 
interval.  The  second  algorithm  takes  these  sets  of  parameters 
and  simulates  the  echo  signal  received  at  the  transducer. 

These  two  algorithms  are  implemented  in  a  single  computer 
program  for  both  M-mode  and  three-dimensional  sector-scanned 
snapshot  simulations. 

Figure  4.1  shows  the  algorithm  used  to  generate  the 
set  of  specular,  scatterer  and  attenuation  parameters.  The 
set  of  geometric  descriptions  is  integral  with  the  program 
for  simplicity  although  the  description  can  be  entered  for 
each  transmit/receive  interval.  First,  the  biological  para¬ 
meters  such  as  velocity  of  sound,  impedance,  attenuation, 
and  surface  scattering  parameters  are  read  in.  For  each 
scan  angle,  or  transmit/receive  interval,  the  boundaries 
are  computed  where  the  cardiac  structure  changes.  At  the 
beginning  of  the  receive  interval,  the  specular,  scatterer 
and  attenuation  parameters  are  set  for  the  whole  interval. 

This  assumes  that  all  subsequent  biological  materials  are 
the  same;  then,  the  range  counter  is  incremented  until  a 
boundary  is  encountered.  At  this  point,  the  specular 
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component  is  computed  from  the  angle  of  incidence  and  the 
characteristic  impedances.  Since  a  relatively  narrow 
beam  of  acoustic  wave  is  used  in  echocardiography,  the 
specular  component  is  further  multiplied  by  the  beam 
function  to  account  for  the  amount  of  energy  reflected  along 
the  beam.  The  backscatter  component  of  the  acoustic  energy 
depends  on  the  roughness  of  the  surface;  this  is  calculated 
according  to  the  characteristic  impedances  and  a  predetermined 
parameter  to  indicate  the  degree  of  roughness.  Since  the 
transmitted  wave  is  attenuated  by  the  boundary,  this 
contribution  is  added  to  the  attenuation  factor  of  the 
particular  biological  material  for  that  range  interval.  The 
attenuation  parameter  used  in  the  signal  simulation  algorithm 
reflects  two  way  propagation  therefore,  the  attenuation 
parameter  set  for  all  subsequent  ranges  is  double  the  incre¬ 
mental  attenuation  factor  for  the  material.  This  process  of 
setting  the  specular  component  at  the  boundary,  and  the 
scatterer  and  the  attenuation  parameters  for  subsequent 
ranges  is  repeated  until  the  end  of  the  receive  interval. 

This  is  repeated  for  the  next  ping  in  M-mode  or  sector-scanned 
mode  until  the  time  or  volume  of  interest  is  completed. 

With  this  set  of  specular  component,  scatterer 
component  and  attenuation  parameters,  the  program  can  generate 
echo  signal  data  according  to  the  signal  model  proposed  in 
Section  2.3.  Figure  4.2  illustrates  the  algorithm  used  to 
generate  echo  signal  data.  The  algorithm  is  a  straightforward 
implementation  of  the  signal  model.  At  each  sampled  range, 
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two  random  numbers  are  generated  for  the  quadrature  compon¬ 
ents.  Both  are  multiplied  by  the  scatterer  component  to 
reflect  the  variance  at  that  range.  The  specular  component 
is  then  added  to  the  cosine  component.  The  total  attenua¬ 
tion  is  updated  and  applied  to  the  quadrature  components. 
Then,  these  components  are  passed  through  different  filters 
to  reflect  the  pulse  shape  and  the  band-pass  process  of  the 
receiver  chain.  Finally,  the  two  components  are  squared, 
summed  and  rooted  to  simulate  the  detector.  Although  the 
actual  acoustic  wave  propagation  and  receiver  processes 
are  more  complicated  than  this  algorithm,  the  power  profile 
is  similar.  This  is  important  in  the  evaluation  of  the 
dynamic  range  compression  algorithm,  and  this  simulation  is 
adequate  for  this  purpose. 

4.1.2  Simulation  of  M-Mode  Echocardiographic 
Data 

The  above  simulation  algorithm  determines 
the  echo  signal  from  a  ping-to-ping  basis.  Therefore,  it  is 
applicable  for  the  simulation  of  an  echo  signal  for  different 
scanning  techniques.  The  only  difference  in  using  the 
simulation  algorithm  for  different  scanning  techniques  is  the 
spatial  or  temporal  sampling  rate,  or,  in  other  words,  the 
geometric  description. 

In  M-mode  echocardiographic  imaging,  the  cardiac 
structural  velocity  encountered  in  an  echocardiographic  scan 
can  be  as  high  as  250  mm/sec.  Therefore,  for  the  typical 
pulse  repetition  frequency  (PRF)  of  1  KHz,  the  simulation 
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model  should  Include  structural  movements  of  about  0.25  mm- 
from  ping  to  ping.  Since  this  rapid  motion  is  generally 
caused  by  aortic  or  mitral  valves,  an  isolated  specular 
reflector  in  a  blood  volume  having  at  least  this  velocity 
is  chosen  for  the  simulation  model.  Figure  4.3  illustrates 
the  structure  used  for  this  simulation  of  M-mode  echocardio- 
graphic  signal  with  dynamic  targets.  The  50  mSec  burst  of 
M-mode  data  contains  velocities  more  than  250  mm/sec,  which 
is  adequate  to  evaluate  the  dynamic  range  compressions 
behavior . 

To  simulate  M-mode  cardiography  with  stationary 
targets,  the  structure  in  Figure  4.4  is  used.  This  is  just 
the  special  case  where  the  targets  are  all  normal  the  acoustic 
beam  in  the  previous  dynamic  M-mode  model.  Although,  the 
targets  are  stationary,  which  implies  the  specular,  scatter 
and  attenuation  parameters  are  the  same  for  different 
pings,  the  actual  echo  signals  for  each  ping  can  differ. 

This  is  due  to  the  nature  of  the  random  process  which  the 
algorithm  simulates,  and  this  reflects  the  actual  signal 
received  under  similar  conditions. 

4.1.3  Simulation  of  Sector-Scanned  Snapshot 
Echocardiographic  Data 

For  a  three-dimensional  snapshot,  it  is  important 
to  include  reflectors  at  various  angles  of  incidence.  The 
set  of  concentric  spheres  shown  in  Figure  4.5  accomplishes 
this  requirement.  It  simulates  the  sector  scanner  as  the 
look  angle  varies.  A  thin  layer  of  fluid  simulates  the 


I,m,:Z,3CT  «  TISSUE 


GO  U. 

N  N 


0  » 


Figure  4.4  Schematic  of  Stationary  Targets  for  M-Mode 
Simulation 
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pericardium.  The  spherical  thin  tissue  within  the  blood 
volume  simulates  valves,  and  the  inner  sphere  of  tissue 
represents  the  intraventricular  septum.  With  this  model, 
the  various  cardiac  structures  are  well  represented. 

4.1.4  Results 

Selected  echo  signal  results  at  various  time 
instants  and  at  various  angles  for  three-dimensional  sector- 
scanned  mode  are  presented  along  with  the  dynamic  range 
compressed  image  in  subsequent  sections.  Table  4.1  summarizes 
the  parameters  used  for  the  biological  materials  for  all  the 
simulations . 

Figure  4.6  shows  the  log  magnitude  of  the  echo 
signal  for  the  stationary  target  case.  Note  that 
particular  ping,  there  is  approximately  a  75  dB  dynamic 
range.  Here,  the  dynamic  range  is  considered  to  be  the  ratio 
of  the  highest  signal  amplitude  of  interest  to  the  signal 
level  of  return  from  blood.  If  the  blood  signal  returns  are 
of  interest,  the  dynamic  range  increases  another  10  dB  or  so. 
At  different  scanning  angles  or  times,  the  blood  signal  level 
might  be  lower  due  to  higher  total  attenuation  in  the  region 
before  the  blood.  Therefore,  the  total  dynamic  range  is 
consistent  with  reported  values  (about  90  to  100  dB) . 

Figure  4.7  shows  the  same  data  as  in  Figure  4.6  in 
linear  scale.  The  first  2.5  cm  is  not  shown  due  to  the  higher 
magnitudes.  Note  that  the  scale  is  arbitrary.  Figure  4.8 
shows  an  actual  A-mode  echocardiographic  signal.  The  low 
signal  level  in  the  first  2.5  cm  is  due  to  the  "long" 
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recovery  time  of  that  particular  amplifier.  Note  the 
similar  nature  of  the  data  after  this  recovery  period  with 
the  simulation  in  Figure  4.6.  The  isolated  specular 
reflectors  and  the  back-scattered  signals  from  tissues  have 
very  similar  appearances  in  the  simulated  data. 

Although  a  direct  prediction  of  the  echo  signal  is 
attempted  here,  it  is  clear  that  simulated  data  are  similar 
in  the  power  profile  to  actual  echocardiographic  data. 

Thus,  the  simulated  data  forms  a  good  basis  for  the  test  of 
the  dynamic  range  compression  algorithm. 

4 . 2  M-Mode  Dynamic  Range  Compression 

A  low-cost  implementation  method  for  the  dynamic 
range  compression  algorithm  is  presented.  The  algorithm  is 
simulated  and  evaluated.  Various  aspects  of  the  algorithm 
are  examined.  A  comparison  with  manual  techniques  of  gain 
control  is  presented. 

4.2.1  M-Mode  Implementation 

In  implementing  an  M-mode  echocardiographic 
imaging  system.  Digital  processing  is  not  most 
economical;  hence,  the  use  of  digitized  echo  data  to 
derive  the  gain  parameter  is  unfeasible.  The  dual  channel 
configuration  suggested  in  Section  3.2.7  is  ideal  for  this 
application.  The  present  analog  video  circuits  require  no 
modification.  The  time  gain  compensation  (TGC)  circuit, 
often  a  separate  module,  can  be  replaced  with  the 
adaptive  gain  unit.  A  schematic  for  such  an  adaptive  gain 
control  is  shown  in  Figure  4.9. 
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The  adaptive  gain  circuit  uses  an  analog  integrator 
to  reduce  the  requirements  on  the  gain  channel  data  rates. 
Instead  of  a  root  mean  square  estimate  of  the  signal  power, 
an  average  estimate  is  used.  By  resetting  the  integrator 
at  the  beginning  of  each  gain  segment  and  using  a  sample  and 
hold  circuit,  the  conversion  rate  for  the  signal  power 
estimate  is  only  12.5  psec/sample.  This  can  be  accomplished 
with  low  cost  analog-to-digital  converters.  Furthermore, 
the  use  of  a  logarithmic  converter  ensures  that  the  dynamic 
range  is  within  the  capability  of  the  analog-to-digital 
converter,  for  example,  48  dB  for  eight-bits.  Since  the 
PRF  for  M-mode  cardiography  is  at  1  KHz  and  the  receive 
interval  for  each  ping  is  about  200  psec . ,  800  p sec.  is 
available  for  deriving  the  gain  needed.  For  simplicity,  a 
sixteen  segment,  constant  gain  profile  can  be  used.  Thus, 
in  between  the  end  of  the  receive  interval  and  the  next 
transmit  times,  the  gain  profile  and  the  bias  parameters  for 
the  next  ping  can  be  derived  and  transmitted  to  the  waveform 
generators.  Figure  4.10  shows  an  example  of  the  timing 
involved. 

The  calculations  involved  in  deriving  the  new  gain 
are  extremely  simple.  First,  in  order  to  utilize  the  ping- 
to-ping  structural  geometric  correlation,  a  simple  digital 
filter  can  be  used.  Then,  for  segment  i,  i-1,.,.,16,  the 
new  power  statistic  can  be  derived  from  the  new  average 
sample  by 


ix. 


A  iX 


n-1 


+  B  iU 


Figure  4.10  Timing  Example  lor  M-Mode  Adaptive  Dynamic 
Range  Compression 
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where  iXn  -  average  estimate  for  segment  i  at  ping  n, 
iUR  *  average  value  for  segment  i  at  ping  n,  and 
A  and  B  are  constants  chosen  for  appropriate  correlation 
time . 

A  minimum  value  is  chosen  to  represent  the 
blood  volume,  min  £iXn^  .  This,  along  with  the  segment 
where  the  minimum  occurs,  can  be  used  to  derive  the 
maximum  and  minimum  bias  levels: 


=  k[  (i  .  -.  5)/5]  2  min 
max  min 


M 


and 


^min 


K  min 
i 


where  K  =  a  constant  depending  on  the  tolerable  false 

alarm  level,  and  i  .  =  segment  where  min  iX  occurs. 

min  n 

It  is  assumed  that  the  near-filed  to  far-field  transi¬ 
tion  zone  occurs  in  the  fifth  segment.  These  two 
parameters  can  then  be  transmitted  to  waveform 
generators  to  form  the  suppression  levels  for  the  next 
ping. 

To  calculate  the  new  gain  parameter  a  multiplier 


is  first  formed: 
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1C  -  max  (kb,  Kd-»  (  min  " - »]  , 

i  1An  7 

where  iC  is  the  multiplier  for  segment  i,  and  where 
K^,  KQ,  and  m  are  determined  from  target  miss  rates 
and  brightness  level  for  isolated  targets  (see  Section 
4.2.5).  Then,  the  new  gain  for  segment  i  can  be 
calculated : 

ig  a  =  S  ViC-iX  , 

6n+l  sat  n  ' 

where 

ign+1  =  new  gain,  and 

Ssat  =  siSnal  saturation  level. 

After  normalization  and  suppression  by  the  gain 
and  bias  parameters,  the  echo  signal  undergoes  a  nonlinear 
transformation  as  expressed  in  Section  4.2.6  to  equalize 
the  distribution.  Specifically,  the  transformation  is: 

T(.)  -  1  -Sjj!  , 

where  N  =  16  for  four-bit  quantization. 


Although  only  5 a  p Sec  are  available  to  derive  the 
new  gain  and  bias  parameters,  careful  design  of  the  hardware 
can  overcome  this  problem.  For  example,  BiUn  can  be  per¬ 
formed  by  setting  the  correct  gain  in  the  analog  amplifiers; 
and  the  division  involved  in  setting  the  new  gain  can  be 
avoided  by  coding  the  waveform  generator  accordingly. 

In  summary,  this  particular  M-mode  implementation 
can  easily  be  implemented  by  using  low-cost  microprocessor 
components.  Figure  4.11  shows  an  off-line  computer 
simulation  of  the  purposed  M-mode  adaptive  dynamic  range 
compression  circuits  used  to  evaluate  performance. 

4.2.2  Results 

The  dynamic -range-compressed  echo -data  for 
stationary  targets  and  dynamic  targets  from  the  M-mode 
simulations  are  presented  here.  The  raw  echo  signal  data 
are  also  presented  for  comparison. 

Figure  4.12  shows  the  log  magnitude  raw  echo 
signal  data  versus  range  from  the  stationary  target  case. 

Note  that  the  raw  signal  amplitudes  are  arbitrary.  The 
absolute  amplitude  depends  on  the  transmitted  power  and  the 
gain  of  the  transducer.  Also,  shown  on  the  same  graph  are 
the  saturation  and  the  bias  (suppression)  levels.  As  seen 
from  the  graph,  both  the  saturation  and  bias  levels  track 
the  signal  fairly  well.  The  isolated  specular  reflectors 
are  shown  to  be  quite  saturated.  But  this  causes  no 
apparent  degradation  of  the  image  because  they  should  be 
displayed  at  the  brightest  levels.  Figure  4.13  shows  the 
intensity  versus  range  of  the  dynamic  range  compressed  image. 
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Figure  4.11  Computer  Implementation  of  M-Mode 
Adaptive  Dynamic  Range  Compression 
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Also  shown  are  the  locations  of  the  different  biological 
materials.  Note  that  the  brightness  level  in  the  blood 
region  is  suppressed.  The  boundary  of  the  different  bio¬ 
logical  materials  can  easily  be  identified.  The  boundaries 
near  the  maximum  range  have  reduced  brightness  levels.  This 
is  due  to  the  lower  limit  on  the  bias  level.  Although  the 
brightness  is  reduced,  the  signal  still  gets  through.  Also, 
since  this  represents  the  far  side  of  the  heart,  performance 
of  the  algorithm  is  not  severely  degraded  due  to  the  relative 
unimportance  in  cardiac  diagnosis.  If  this  region  is  of 
special  interest,  such  as  the  detection  of  pericardial 
effusion,  the  operator  can  set  the  bias  level  lower  just  for 
this  particular  operation  only.  In  general,  the  algorithm 
produces  the  image  desired. 

Figure  4.14  shows  the  schematic  of  the  M-mode  dynamic 
target  image  simulated.  Three  pings  are  selected  from  this 
set  to  show  the  result  of  the  algorithm.  Figures  4.15  to 
4.17  show  the  log  magnitude  of  the  raw  data  and  the  intensity 
versus  range  data  for  the  three  examples.  The  component  for 
the  stationary  target  applies  to  this  image  data  also. 

4.2.3  Discussion 

Several  performance  aspects  of  the  M-mode  implementation 
of  the  dynamic-range-compression  algorithm  are  presented  here. 

In  Section  4.2.2,  we  see  that  the  algorithm  yields  the  image 
we  desired.  In  this  section,  aspects  such  as  intensity  dis¬ 
tribution,  transient  behavior,  effects  of  quantization  on 
power-statistic  estimation  and  gain  commands  as  well  as 
compression  with  manual  gain  control  are  discussed. 
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Figure  4.14  Schematic  of  M-Mode  Moving  Targets 
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Figure  4.15  M-Mode,  Moving  Targets,  Example  1 
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4.2.3. 1  Symbol  Distribution 

Ideally,  the  distribution  of  symbols  should  be 
uniform  to  maximize  the  information  and  to  enhance  the 
echocardiographic  image.  There  are  factors  which  will 
cause  the  distribution  to  deviate  from  uniform.  First, 
the  null  symbol,  used  to  represent  blood,  is  dependent  on 
the  amount  of  blood  volume  for  that  particular  scan.  In 
other  words,  the  percentage  of  the  symbol  zero  in  a  scan 
is  image  dependent.  The  same  applies  for  the  maximum 
symbol.  If  enough  isolated  specular  reflection  occurs  in  a 
scan,  then  there  will  be  proportionately  more  signals  of 
maximum  amplitude.  In  order  to  investigate  the  symbol  dis¬ 
tribution,  the  bias  parameter  should  be  neglected  and  the 
scan  should  be  relatively  free  of  isolated  specular  reflectors. 

The  other  factors  that  affect  the  symbol  distribution 
are  due  to  estimation  errors  and  the  implementation  method. 
Since  the  gain  parameters  are  applied  for  a  small  segment, 
namely,  one-sixteenth  of  the  receive  interval,  then  due  to 
the  nonstationarity  of  the  signal,  an  estimation  error  will 
occur  on  the  signal  statistics.  Also,  since  all  points 
within  this  segment  are  applied  with  the  same  gain  parameter, 
all  points  cannot  be  normalized  correctly.  Therefore,  the 
symbol  distribution  will  deviate  from  the  ideal  uniform 
distribution . 

Figure  4.18  shows  a  typical  M-mode  scan  with  the  gain 
parameter  but  not  the  bias  parameter  of  the  dynamic  range 
compression  algorithm  applied.  This  ping  contains  a  region 
of  blood  and  tissue  but  does  not  contain  isolated  specular 
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Figure  4.18  Distribution  of  Symbols  (Intensity) 
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reflectors.  The  top  graph  shows  the  log  magnitude  signal 
data  as  a  function  of  range.  Also,  the  saturation  level 
is  shown.  Since  the  gain  profile  is  derived  from  the  past 
pings  in  the  M-mode  implementation,  the  gains  derived  are 
not  perfect,  but  obviously  adequate.  The  middle  graph 
shows  the  echo  signal  after  adaptive  gain  and  quantization. 
Note  that  the  signal  appears  homogeneous  and  noisy  throughout 
the  scanning  range.  This  illustrates  the  importance  of  the 
bias  parameter.  In  order  to  automatically  apply  gain  to  the 
echocardiographic  image,  the  suppression  level  must  also  be 
coupled  so  that  the  low  level  signal  from  a  blood  volume 
will  not  be  automatically  magnified  and  thus  render  the  image 
useless.  Also,  note  that  since  the  segment  transition 
points  do  not  occur  exactly  at  the  signal  transition  points, 
the  high  and  low  magnitude  regions  are  separated  by  a  dark 
band.  This  results  in  additional  information  to  the  operator 
to  signify  a  power  transition  point.  This  is  similar  to  edge 
enhancement  by  differentiation.  The  bottom  graph  shows  the 
resultant  histogram.  Even  though  the  gain  profile  is  derived 
from  past  pings,  the  histograms  are  still  fairly  uniform  due 
to  the  high  correlation  from  ping  to  ping.  Clearly,  a  fairly 
uniform  symbol  distribution  can  be  achieved. 

4. 2. 3. 2  Transient  Behavior 

In  order  to  increase  the  estimation  accuracy  in  the 
power  statistics,  a  digital  filter  is  used  from  ping  to  ping 
for  the  same  segments  in  the  dynamic  range  compression 
algorithm.  Figure  4.19  shows  the  transient  response  of  the 
filter  used  in  this  particular  implementation  of  the  dynamic 
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range  compression  algorithm.  Although  the  correlation  time 
of  the  cardiac  structure  is  longer,  the  shorter  transient 
response  speeds  up  the  response  to  transient  conditions. 

For  example,  with  structure  velocity  of  0.25  mm/ping,  then 
it  will  take  3.75  pings  for  the  structure  to  move  across  the 
segment.  Since  there  are  essentially  25  independent  samples 
available  per  segment  for  the  power  statistic  estimate,  it  is 
sufficient  to  use  a  faster  filter  for  better  transient 
response. 

Figure  4.20  shows  how  the  sixteen  gain  segments 
converge  to  the  final  valve  within  two  pings.  This  is  for 
the  M-mode  simulation  of  a  stationary  target.  Note  that  the 
variation  after  the  third  ping  is  quite  small.  Some  varia¬ 
tion  can  be  attributed  to  the  fluctuation  of  the  signal  it¬ 
self;  a  stationary  target  means  constant  specular,  scatterer 
and  attenuation  parameters  only,  but  not  a  constant  signal. 

From  the  convergence  rate  of  the  algorithm  shown,  it 
is  clear  that  it  can  handle  rapid  signal  variations  from  ping 
to  ping.  This  digital  implementation  for  frame  to  frame 
operation  has  clear  advantages  over  analog  techniques.  If 
the  operator  wishes,  to  change  the  filter  constant  in  this 
digital  implementation,  it  involves  changing  three  memory 
words  only.  Therefore,  to  utilize  the  inherent  signal 
correlation  from  frame  to  frame,  the  digital  implementation 
is  adequate  to  derive  an  adaptive  gain  for  echocardiographic 
imaging . 
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4. 2. 3. 3  Comparison  with  Manual  Gain  Control 

Because  the  manual  gain  control  used  in  present 
day  echocardiographic  imaging  systems  depends  on  operator 
control,  an  objective  direct  comparison  with  the  adaptive 
dynamic  range  compression  algorithm  is  impossible.  However, 
one  can  compare  the  adaptability  of  the  different  gain 
control  systems  with  dynamic  targets. 

A  percentage  error  approach  is  used  here  for  several 
reasons.  Fidelity  criteria  for  images  are  hard  to  define 
objectively  and  consistently  for  different  types  of  images. 
Furthermore,  a  transformation  from  raw  echo  data  to  image 
data  cannot  be  readily  defined.  The  percentage  error 
approach  applies  if  one  wishes  to  reconstruct  the  full 
dynamic  range  echo  data  by  applying  the  inverse  gain.  This 
approach  has  applications  if  the  full  dynamic  range  data  is 
used  in  further  processing  or  enhancement,  for  example,  in 
beam  width  deconvolution  or  compound  scanned  image  reconstruc¬ 
tion.  Figure  4.21  illustrates  the  process  used  to  evaluate 
the  average  percentage  error  for  the  various  gain  control 
processes. 

The  manual  gain  control  methods  used  in  existing 
echocardiographic  systems  are  constant  throughout  the  image. 
Therefore,  it  can  only  be  optimized  for  a  single  image  or 
ping.  The  actual  gain  applied  depends  on  the  operator  skill. 
Therefore,  the  percentage  error  presented  here  only  serves  as 
comparison  between  stationary  target  and  the  effects  of  the 
target  motion.  For  the  gain  control,  the  parameters  are 
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derived  from  a  stationary  target,  consistant  with  operating 
procedure.  Figure  4.6  illustrates  the  echo  signal  selected 
for  this  purpose.  The  log  magnitude  data  is  used  to  derive 
the  gain  parameter.  Usually,  a  static  log  magnitude  display 
is  not  available  to  the  operator;  thus,  this  serves  as  a 
better  guide  in  deriving  the  gain  controls. 

A  logarithmic  compression,  a  single  slope  time  gain 
compensation  (TGC)  and  slope  and  logarithmic  manual  gain 
control  methods  are  selected  as  the  baseline  to  compare  to 
the  adaptability  of  the  dynamic  gain  control  algorithm. 

Their  corresponding  schematics  and  the  gain  control  parameters 
are  presented  in  Figure  4.22.  A  summary  of  the  average  per¬ 
centage  errors  for  the  three  methods  for  the  constant  target 
and  the  dynamic  target  cases  are  presented  in  Table  4.2. 

Since  reconstruction  of  the  full  dynamic  range  data 
is  desired,  the  bias  parameter  in  the  dynamic  range  compression 
algorithm  is  set  to  zero  for  this  purpose.  Figure  4.23  shows 
the  average  percentage  errors  for  the  adaptive  dynamic  range 
compression  algorithm  for  various  values  of  the  parameter  K^. 
Note  that  the  errors  reach  a  minimum  for  some  value  of  K 

00 

and  is  relatively  insensitive  around  the  minimum  value.  The 
shape  of  the  curve  is  expected  because  the  error  consists  of 
three  terms:  quantization  error,  saturation  error  and  low 
level  exclusion  error.  Although  the  latter  two  are  mani¬ 
festations  of  quantization,  they  are  highly  dependent  on  the 
probability  density  of  the  signal.  Since  the  low  level  error 
has  a  percentage  error  of  100%  and  the  saturation  error  is 
fairly  high,  we  can  assume  the  form  of  the  total  error  to  be 


LOG  COMPRESSION: 


y  =  4^  LOG  108x  ;  ye  [0. 1  ] 


SINGLE  SLOPE  GAIN: 


r  =  RANGE  (in cm) 


SINGLE  SLOPE  GAIN  WITH  LOG  COMPRESSION: 


Figure  4.22 
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L  ~U 

STotal  “  EQ  +  j  PS(s)ds  +  j  PS(s)ds  ’ 


where 


Eq  =  quantization  error, 

S^  =  low  level  exclusion  level  ,  and 
Sjj  =  saturation  level. 


SL  and  Sy  are  functions  of  the  number  of  quantization  levels 
and  of  the  inverse  nonlinear  transformation  T_1(s) .  Since 
the  gain  applied  is  assuming  aQ=o  ,  we  have  ; 

2  2 
2  “Cs 

PS(S)  =  KooS  ® 


and  therefore, 


ETotal  EQ  +  ^  1-6 


-st2k2 
L  x  j 


2  2 
-S  ,K* 
u  x 


"SLKx 

The  expression  (1-e  )  increases  as  K  increases,  while 

2  2  X 

ux 

e  decreases  as  K  increases.  Since  S*»S.T,  the  form  of 

x  L  U 

the  total  error  as  a  function  of  is  concave  as  in  Figure 

4.23.  For  between  1.75  and  2.75,  the  total  percentage 

error  changes  very  little.  This  insensitivity  helps  the 

designer  in  choosing  K  . 

x 

For  =  2.25,  the  total  percentage  error  for  the 
adaptive  dynamic  range  compression  algorithm  for  stationary 
targets  and  for  moving  targets  are  presented  in  Table  4.2. 


Table  4.2  Percentage  Errors  for  Different 

Dynamic  Range  Compression  Algorithm 


STATIONARY 

TARGET 

MOVING 

TARGET 

% 

CHANGE 

LOG  COMPRESSION 

26.6  dB 

36.5  dB 

37.2  dB 

SINGLE  SLOPE 

23.7 

34.5 

45.6 

SLOPE  W/LOG 

16.4 

25.7 

56.7 

ADAPTIVE  COMP. 


17.6 


20.4 


15.9 
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Note  that  the  single  slope  with  logarithmic  compression  has 
the  lowest  percentage  error  for  stationary  targets.  This  is 
expected  because  logarithmic  compression  yields  the  lowest 
percentage  error  no  matter  what  the  distribution  is.  Since 
the  dynamic  range  compression  algorithm  is  designed  for 
imaging,  the  logarithm  compression  approach  was  not  used. 

The  percentage  change  from  stationary  targets  to 
moving  targets  for  the  various  gain  controls  are  also  shown 
in  Table  4.2.  Obviously,  the  dynamic  range  compression  pre¬ 
sented  has  the  best  performance  in  a  dynamic  target  environ¬ 
ment.  Also,  note  that  no  operator  adjustment  is  required 
even  for  one  single  scan. 

4. 2. 3. 4  Effects  of  Quantization 

An  actual  implementation  of  the  M-mode  adaptive 
dynamic  range  compression  algorithms  will  utilize  an  eight-bit 
A/D  converter  for  the  purpose  of  estimating  the  power  statis¬ 
tics  and  will  utilize  an  eight-bit  control  word  to  the  voltage 
control  amplifier.  An  examination  of  what  effect  this  will 
have  on  the  image  will  be  examined  here.  Again,  the  average 
percentage  error  is  used  to  evaluate  this  effect.  Theoretic¬ 
ally,  the  quantization  effects  are  minimal.  Quantization 
adds  a  random  error  term.  For  eight-bit  quantization,  this 
introduces  a  noise  level  58.9  dB  below  the  signal  level. 

Even  if  this  error  is  fairly  high,  the  relative  insensitivity 
of  the  algorithm  to  Kx  will  not  degrade  the  saturation  per¬ 
formance  much.  The  suppression  level  will  suffer  slightly 
due  to  quantization.  For  eight-bit  quantization,  this  degrada¬ 
tion  is  negligible.  Figure  4.24  shows  the  average  percentage 
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error  with  and  without  quantization  in  the  power  estimate 
and  gain  control  versus  K^.  Obviously,  the  effect  is 
negligible. 

4 . 3  Three-Dimensional  Sector-Scanned  Snapshot 

Mode  Dynamic  Range  Compression 

For  multidimensional  echocardiography,  the  dynamic 
range  compression  algorithm  may  be  implemented  differently  from 
the  M-mode  technique  suggested.  One  reason  is  the  inherent 
complexity  of  multidimensional  scanners.  The  addition  of 
hardware  or  software  to  implement  the  dynamic  range  compression 
algorithm  is  of  little  impact  to  the  system.  Such  an  example 
is  the  three-dimensional  sector-scanned  snapshot  camera 
suggested  by  Miller,  Lawther  and  Zelis.^9^  The  snapshot 
camera  utilizes  digital  memory  for  the  bulk  storage  of  three- 
dimensional  sector-scanned  echocardiographic  signal  data. 
Subsequent  reconstruction  of  sector  display  is  controlled  by 
a  minicomputer.  The  addition  of  control  interfaces  for  volt¬ 
age  control  amplifiers  and  digital  look-up  tables  do  not 
increase  the  system  complexity  by  a  significant  amount. 

The  following  snapshot  mode  dynamic  range  compression  algorithm 
is  designed  specifically  for  this  echocardiographic  system. 

4.3.1  Snapshot  Mode  Implementation 

For  three-dimensional  snapshot,  the  dynamic  range 
compression  algorithm  can  be  applied  after  data  acquisition. 
Since  an  eight-bit  ana log-to-digital  converter  is  used, 
logarithmic  compression  and  subsequent  expansion  should  yield 
reasonable  data  with  post  scan  expansion.  With  the  full 
reconstructed  dynamic  range  data,  one  can  apply  the  dynamic 
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range  compression  algorithm  as  suggested  in  Chapter  III. 

For  a  large  amount  of  data,  the  post  scan  data  processing 
time  can  be  quite  long,  therefore,  an  alternative 
application  of  the  dynamic  range  compression  algorithm  with 
partial  hardware  implementation  is  desirable.  Figure  4.25 
shows  a  schematic  of  such  a  design. 

The  operational  aspects  of  the  snapshot  camera  allows 
the  possibility  for  an  initial  scan  to  acquire  the  gain  and 
bias  parameters.  Then,  the  actual  diagnostic  scan  can  be 
taken.  This  single  channel  realization  eliminates  the  dual 
channel  electronics.  Since  a  high-speed  analog-to-digital 
converter  is  used  in  conjunction  with  a  digital  memory  for 
the  snapshot  camera,  the  dual  channel  slow-conversion 
approach  for  the  M-mode  imaging  systems  seems  unnecessary. 

In  Figure  4.25,  the  receiver  consists  of  a  voltage 
controlled  amplifier  detector,  logarithmic  amplifier,  eight- 
bit  A/D  converter,  a  programmable  table,  and  waveform  genera¬ 
tor.  The  digital  memory  provides  data  storage,  and  the  CPU 
is  used  to  calculate  gain  and  bias.  The  CPU  also  controls  a 
grey  scale  display  unit  which  uses  the  data  from  memory. 

The  waveform  generator  provides  a  piece-wise  linear  continuous 
gain  profile  (in  dB)  to  the  UCA .  Sixteen  instruction  words 
provide  the  initial  gain  and  incremental  gains  as  illustrated 
in  Figure  4.26.  The  piece-wise  linear  continuous  gain  profile 
eliminates  the  abrupt  switching  transient  as  in  the  M-mode 
implementation.  Also,  the  manually  adjusted  piece-wise 
linear  continuous  gain  profile  is  used  clinically;  thus, 
acceptability  of  the  gain  profile  by  the  end  user  is  enhanced. 


Figure  4.25  3D  Snapshot  Implementation  of  Adaptive  Dynamic  Range 
Compression  Algorithm 
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The  voltage-controlled  gain  amplifier  should  be  located 
at  the  front  end  of  the  receiver  chain  to  keep  the  signal- 
to-noise  ratio  of  the  latter  stages  high.  Since  the  bias 
parameter  is  not  critical,  the  bias  parameter  can  be 
applied  digitally  and  in  real  time  with  a  table  look-up 
technique.  A  programmable  table  is  connected  to  the  eight- 
bit  A/D  converter  and  updated  as  the  range  changes.  This 
will  dynamically  apply  the  bias  versus  range  profile.  The 
effects  of  segmenting  the  bias  parameter  are  minimal.  If 
sixteen  segments  are  used,  the  first  segment  is  about  6  cm 
long  in  range  for  the  near  zone.  The  remaining  segments 
are  0.6  cm.  Therefore,  the  largest  error  in  bias  will 
occur  adjacent  to  the  transition  zone.  This  will  contribute 
10%  error  to  the  bias  parameter  at  this  maximum  error  range. 
Since  the  bias  is  a  lower  limit  for  blood  level  only,  this 
10%  error  will  not  affect  the  suppression  function.  Of 
course,  at  higher  ranges,  this  error  is  smaller,  the  minimum 
being  at  the  maximum  range  with  4.3%.  This  programmable 
table  also  provides  the  necessary  nonlinear  transformation 
to  the  log-compressed  data  to  yield  a  uniform  output  distri¬ 
bution.  Since  the  transformation  is  programmable,  one  can 
use  a  different  one  as  desired.  This  configuration  offers 
tremendous  flexibility.  A  detailed  schematic  of  this 
programmable  table  is  shown  in  Figure  4.27.  The  main  function 
of  the  table  is  to  form  a  four-bit  word  from  an  eight-bit 
input.  This  can  be  accomplished  by  using  random  access 
memories  (RAMs) .  The  eight-bit  input  forms  part  of  the 
address  to  the  RAMs.  Another  four  bits  forms  the  address  to 


Figure  4.27  Look-up  Table  Schematic 
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advance  the  table  as  a  function  of  range.  A  clock  and 
counter  can  be  set  to  increment  range  counts  at  the  far 
zone.  The  four-bit  output  words  can  then  be  packed  into 
eight  bits  so  that  it  is  compatible  with  the  interface. 

At  each  of  the  sixteen  bias  range  increments,  the  bias 
level  can  occur  at  any  one  of  the  lower  levels.  Figure  4.28 
illustrates  the  bias  level  for  different  segments.  For  this 
illustration, 64  levels  are  suppressed  by  the  bias  level; 
there  are  192  levels  left  over  to  perform  the  transformation. 
Figure  4.29  illustrates  the  eight-bit  input,  four-bit  output 
transformation.  Under  actual  operating  conditions,  since  the 
bias  parameter  and  the  gain  parameter  are  uncorrelated  for 
nonblood  volumes  and  correlated  for  blood  volumes,  the  bias 
parameters  can  increase  or  decrease  from  one  segment  to 
another.  Therefore,  the  remaining  range  encoding  into 
nonzero  symbols  also  changes.  This  necessitates  a  new  trans¬ 
formation  for  that  particular  segment.  Therefore,  the  six¬ 
teen  look-up  tables  vary  for  each  of  the  bias  segment  for 
each  statistical  window.  This  implementation  of  a  look-up 
table  enables  the  real-time  application  of  the  algorithm  by 
eliminating  the  point-by-point  computations  required. 

The  application  of  the  above  implementation 
technique  consists  of  the  following  steps  and  calculations: 

1)  Transfer  initial  gain  vector  to  buffer. 

Let  G  =  [ Gq , , . . . , G^g] ^  , 
where 

GQ  -  initial  gain  in  dB,  the  selection  of 


m 


Figure  4.29  Continued 
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Figure  4.29  Continued 


which  depends  on  the  approximate  initial 
signal  strength. 

G1,...,G15  =  +2  dB/cm  (for  example) 

=  incremental  gain  in  dB, 

2  dB/segment  is  chosen  to  compensate  for 
some  attenuation  effects  but  low  enough 
so  that  the  latter  segments  are 
saturation  free. 

Commence  gain  parameter  acquisition  scan. 

The  eight-bit  log-compressed  data  for  a 
complete  scan  are  stored  in  memory  for  the 
derivation  of  gain  and  bias  parameters. 

Let  x .  .  .  =  log  g .  s .  .. 

where 

x.  .  .  =  log  compressed  data  as  a  function 

1  >  J  >  K 

of  range  and  scan  angles  (see 
Table  4.2.2. 1) 

gi  =  initial  gain  for  range  i. 
s.  .  ,  =  raw  echo  signal. 

1  >  J  >  K 

Derive  gain  and  bias  parameters  for  diagnostic 
scan.  The  calculations  are  as  follows: 
a)  Calculate  variance  for  each  statistic 
using  correct  window  and  range.  The 
window  refers  to  box  car  averages  over 
the  two  scanning  angles 
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los  si,j,k  ‘  xi,jk 


s?,jk  ‘  l0«'112  lo*  si,J,k] 


-i.p.,  -japr  ZZE  =2  • 

I  j  K  BOX  CAR  SUM 


°i,P,q 


where  p  =  index  of  angle  segment 
over  scanning  plane  9,  and 
q  =  index  of  angle  segment  over 
scanning  plane  0. 


Also  let 


o  .  =  min ja.  ? 

min  i ,  P  i  q  3 


and 


min 


b)  Calculate  gain  segments, 


^i,p,q  f<Kl’V  °i,p,q  °min^ 


where 
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f  =  nonlinear  function  to  limit  signal 
near  the  minimum  level  (see  Section 
3.2.5) . 


gi , p, q 


1/K . 


a. 


i,p,q  i,p,q 


9 


A 

gi, p, q 


=  20  log  g 


i,p,q 


9 


where  g.  =  new  gain  at  range  i, 

1  9  P  9  4 

angle  segments  p,q,  and 

in  dB's. 


=  g 

&i , p, q  i ,  p,  q 


Fit  g  to  piece-wise  linear  continuous  curve  by 


T 

(FAF) 


A 

F  g. 

Jci,  p,q 


» 


where  F  =  f(i)  =  function  to  generate  piece- 
wise  linear  continuous  function,  and 
G  =  new  gain  parameters  for  angle 

~p,q 

segments  p,q. 
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After  all  gain  parameters  are  used  to  solve 
for  the  different  angles  and  angle  segments, 
transmit  to  buffer. 

c)  Calculate  bias  parameter 

First  calculate  average  gain  over  segment, 
(since  bias  refers  to  the  raw  signal  and  in 
this  implementation  the  suppression  is 
applied  post  multiplication,  and  multiplica¬ 
tive  constant  must  be  removed.) 

g  =  ^  WT  G 

sn ,  p,  q  20  ~n  ~p,q 

where  g  =  average  gain  in  dB's  over 

n  f  P  t  Q. 

bias  segment  n,  and  JH  =  weighting 
function  for  bias  segment  n, 


The  minimum  and  the  maximum  bias  parameters  can 
be  calculated  by, 


^min 

^min 


Kb  “min 

/  min  ^  u 
<7^~  "i" 


min 


*min 

> 

*FF 

*min 

£ 

lFT 
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where  Kb  =  constant  dependent  on  missed 
target  rate,  and  IFF  -  far-field  range. 


Then  the  segmented  and  quantized  gain  applied 
with  the  look-up  table  can  be  calculated  by: 


B  -  g  +  log  b 

o,p,q  so,p,q  6  max 


1  2 

B  =  g  +  max  b  .  ,  (— — )  b 

n,p,q  Bn, p,q  min’  IFF  max 


for  n  =  1,2, . . . , 15, 


where  In  =  average  range  for  bias  segment  n. 


The  quantized  value  of  B  (in  this  case 

n,p,q 

[o,  ^»,  .  .  .  ] )  can  then  be  used  to  set 

the  look-up  table  zero  output  symbols. 

d)  Calculate  breakpoints  for  look-up  table 

The  nonzero  ranges,  B  to  1,  should  be 

n  9  P » 4 

grouped  among  fifteen  symbols  according  to  the 
cumulative  distribution  transformation.  Since 


the  cumulative  distribution  transformation  is 


where  N  -  16. 


And  the  log  amplifier  transformation  is: 


y  -  3  log  x+1 


The  breakpoints  are: 


'i  ‘  Wo-°.,p Bn,p,q1+l, 


where  i  =  1,2,..., 15,  and  *  decision 
levels  for  the  look-up  table  transforma 
tion. 


e)  Set  look-up  table 

With  B  and  its  corresponding  y.,  the 

11  9  P  f  4  1 

output  symbols  can  now  be  set  at  decision  points 

^Bn,p,q’  yi»  y2’*'',y15^  with  corresponding 
symbols  (0 , 1 , 2 , . . . , 15) . 
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4)  Commence  diagnostic  scan 

As  scanning  plane  and  angle  increments,  the 
different  gain  and  bias  parameters  are 
applied  by  the  waveform  generator  and  multi¬ 
plier,  and  by  the  RAM  look-up  table, 
respectively . 

As  shown  in  the  above  algorithm,  the  hardware 
takes  over  the  post  processing  requirement  for 
the  computer.  In  an  actual  hardware  implementa¬ 
tion,  the  buffers  can  be  eliminated.  The  same 
computer  memory  used  for  storing  the  initial 
data  and  the  final  image  can  be  used  for 
storing  the  gain  and  bias  parameters.  The  data 
rate  for  the  final  image  is  reduced  to  eight-M- 
bits/sec  assuming  a  2  MHz  sampling  rate  of  the 
ultrasonic  signal.  This  is  well  within  the  one 
M-words/sec  DMA  transfer  rate  of  a  16-bit/word 
minicomputer.  Therefore,  time  can  be  shared 
between  storing  the  image  data  and  transferring 
the  gain  and  bias  parameters  through  the  DMA 
channel  of  the  minicomputer.  Furthermore,  this 
technique  can  easily  adapted  to  real  time 
sector  scanning.  Since  ample  data  samples  are 
available  for  each  statistic  window,  these 
samples  can  be  reduced  in  the  initial  scan  so 
that  data  over  time  can  also  be  sampled  to 
derive  a  set  of  gain  and  bias  parameters  over 
four  dimensions. 


4.3.2  Results 


Figure  4.29  shows  the  computer  algorithm  to  emulate 
the  three-dimensional  sector-scanned  snapshot  mode  adaptive 
dynamic  range  compression  algorithm  as  outlined  in  the  last 
subsection.  The  program  reads  in  a  complete  sector  scanned 
snapshot  during  the  first  pass  and  uses  it  to  derive  the 
gain  parameters  and  the  look-up  tables.  During  this  pass, 
an  initial  gain  profile  and  a  logarithmic  transformation  are 
used  to  compress  the  echo  data  to  eight  bits,  simulating  the 
actual  process  of  the  snapshot  earner  hardware.  The 
statistic  window  used  to  derive  these  parameters  is  three 
pings  by  three  pings,  corresponding  to  a  three  degrees  solid 
angle.  After  the  gain  profile  and  bias  parameters  are 
calculated  for  that  statistic  window,  the  gain  profile  in  dB 
verse  range  is  fitted  to  a  piece-wise  linear  continuous  curve 
and  the  bias  parameters  are  used  to  set  up  the  look-up  table. 
The  program  is  then  ready  for  the  second  pass.  The  input 
device  is  rewound.  As  each  ping  of  the  three-dimensional 
sector  scan  is  read  in  during  the  second  pass,  the  gain  pro¬ 
file  is  applied  and  the  result  is  log-compressed  and 
quantized.  The  eight-bit  log-compressed  signal  is  then  used 
in  the  look-up  table  to  apply  the  bias  and  the  histogram 
equalization  transformations.  The  four-bit  data  is  now  ready 
for  plotting  or  other  evaluation  purposes. 

Figure  4.30  shows  a  typical  set  of  raw  and  compressed 
data.  Also  shown  is  the  saturation  level.  Although  the 
saturation  level  does  not  match  the  raw  data  exactly  for  this 
ping,  it  clearly  is  adequate  for  tracking  the  power  profile. 


From  the  adaptive  dynamic  range  compressed  graphs,  it  is 
clear  that  the  algorithm  distinguishes  blood  volume  with  a 
minimum  of  artifacts  and  at  the  same  time  offers  varied 
intensity  levels  in  arrears  with  significant  returns.  This 
particular  implementation  is  similar  to  an  M-mode  technique; 
therefore,  all  discussion  of  the  algorithm  also  applies  to 
this  technique. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

In  echocardiography,  there  are  two  types  of 
signals  present  at  the  receiver.  One  is  from  specular 
returns  due  to  abrupt  changes  in  acoustic  impedance; 
the  other  is  scatterer  returns  due  to  heterogeneous 
tissue  layers  or  blood  cells.  As  an  echocardiographic 
image,  the  specular  returns  and  the  scatterer  returns  from 
tissues  are  of  interest.  Therefore,  for  a  grey  scale 
display,  maximum  brightness  for  the  specular  reflectors 
and  a  grey  scale  for  the  tissue  scatterers  form  an 
acceptable  image.  In  order  to  generate  an  image  from 
the  wide  dynamic  range  echo  data,  the  image  processor  must 
adaptively  compensate  for  the  attenuation  effects  and 
automatically  reject  returns  from  a  blood  volume.  Vari¬ 
able  gain  and  bias  parameters  are  suggested  to  adaptively 
normalize  the  desired  echo  signals  and  to  adaptively 
reject  the  undesired  echo  signals.  Two  economically 
feasible  implementation  techniques  were  shown.  A  model 
and  its  simulated  echo  data  were  used  to  evaluate  the 
above  techniques  for  M  mode  and  three-dimensional  snapshot 
echocardiography.  Performance  for  both  techniques  was 
found  satisfactory;  no  other  adaptive  techniques  are 
available . 

For  future  research,  many  improvements  can  be  made 
to  the  echo  signal  model  and  adaptive  dynamic  range 
compression  algorithm.  The  model  developed  herein  is 
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primarily  for  the  test  of  the  adaptive  algorithm  response 
to  different  power  profiles.  But  the  model  has  shown 
promise  to  simulate  echo  data  if  exact  cardiac  geometry 
and  biological  parameters  are  used  as  inputs.  In-vitro 
experiments  with  fairly  complicated  structures  can  be 
performed  to  verify  the  echo  signal  model.  This  will 
enhance  an  understanding  of  the  available  single-structure 
biological  parameters.  With  this  verification  of  model, 
perhaps  an  iterative  process  of  regressing  the  actual  echo 
cardiac  signal  data  to  fit  a  cardiac  structure  model  is 
feasible.  This  process  is  in  accordance  with  the  ultimate 
goal  of  identifying  different  cardiac  structures 
automatically.  Another  way  to  achieve  this  goal  is 
through  the  use  of  pattern  recognition  techniques.  With 
the  adaptive  dynamic  range  compression  technique,  full 
dynamic  range  echo  data  can  be  reconstructed  with  off-line 
computers.  Since  blood  volumes  have  lower  amplitude  and 
higher  doppler  shifts  than  other  cardiac  structures  in 
general,  pattern  recognition  techniques  can  be  applied  on 
these  two  features  to  extract  cardiac  chamber  geometries. 
This  two-feature  pattern  recognition  technique  should 
offer  better  performance  than  the  single-feature  (level 
discrimination)  technique  suggested  here. 

In  summary,  the  adaptive  dynamic  range  compression 
algorithm  is  capable  of  compressing  raw  echo  data  directly 
to  display  level  dynamic  ranges.  Also,  since  the  technique 


allows  the  digital  storage  of  echocardiographic  data 
economically  and  permits  the  reconstruction  of  the  full 
dynamic  range  raw  echo  data,  it  opens  endless 
possibilities  to  process  the  echo  data  to  suit  diagnostic 


needs . 
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APPENDIX  A 


HUMAN  TESTS  OF  ADAPTIVE  TVG 
USING  M-MODE  ECHOCARDIOGRAPHY 

In  this  appendix,  one-dimensional  dynamic  range 
compression  of  real  A-mode  echocardiographic  data  is 
presented.  The  data  was  taken  from  normal  volunteers 
with  a  commerical  echocardiographic  imaging  system.  The 
received  echo  signal  was  also  fed  to  manually  adjust 
multisegment  time  varying  gain  (TVG)  controlled  amplifiers. 
This  manually  dynamic  range  compressed  signal  was 
subsequently  recorded  on  digital  magnetic  tape  for  off¬ 
line  analysis.  Since  the  amplifiers  were  calibrated  and 
the  gain  control  signal  were  recorded,  an  off-line 
reconstruction  of  the  full  dynamic  range  echo  signal  was 
possible.  This  full  dynamic  echo  data  were  then  compressed 
with  the  stated  algorithm.  Such  a  technique  can  be  applied 
to  the  three-dimensional  snapshot  camera.  But  due  to  the 
technique  presented  in  Section  4.3.  This  appendix  also 
includes  results  from  the  off-line  implementation  of  the 
dynamic  range  compression  algorithm  operating  of  real  data. 

Medical  diagnostic  pulse-echo  instruments  produce 

(69) 

echo  signals  with  dynamic  range  of  about  100  dB,  ' 
depending  upon  the  application.  Echocardiographic  signals, 
which  typically  span  dynamic  ranges  of  90  dB,  need  to  be 
compressed,  with  minimum  resulting  image  distortion.  The 
requirement  on  compression  of  dynamic  range  is  set  mostly 
by  display  devices,  which  are  limited  to  less  than  15  dB. 
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Several  two  dimensional-echocardiographic 
scanners  and  experimental  results  are  desired  in  the 
literature. (*5,61)  The  transducer  apertures  typically 
vary  from  10  mm  diameter  to  arrays  12  mm  by  24  mm.  The 
reasonant  transducer  frequency  range  is  primarily  from 
2.25  to  2.5  MHz.  Two-  and  three-dimensional  scanners 
are  being  studied,  including  scanning  pistons  and  a  two- 
dimensional  electronically-phased  transducer  array. 

The  basic  elements  of  a  pulse-echo  system  were 
shown  in  Figure  1.2.  Large  diameter  ranges  are  the 
result  of  ultrasonic  absorption,^68)  reflection,  and  wave 
divergence.  Echo  power  profiles,  being  dependent  on 
tissue  structure,  differ  with  transducer  look  direction 
and  with  the  phase  of  the  cardiac  cycle.  These  cannot  be 
compensated  by  operator  manual  adjustment  in  sector¬ 
scanning  echocardiography.  Repeatable  results  cannot  be 
assured  using  manual  adjustments  of  time  varying  gain  (TVG) . 

Also,  a  major  workload  is  involved  in  control  manipula- 
(7) 

tion. v  These  considersations  have  motivated  work  on 

adaptive  TVG  profiles  or  AGC  (adaptive  gain  control)  for 

ultrasound  imaging.  AGC  systems  for  many  signal  processing 

applications  previously  have  been  developed  and 

studied. ^ ^ For  example,  DeClercq  and 
(7 ) 

Maginess'  describe  an  analog  AGC  for  dynamic  ultrasound 
imaging  and  discuss  the  relationship  to  manual  ultrasonic 
TVG  and  radio  AGC;  Stockham  describes  the  application  of 
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(58) 

homomorphic  filtering  to  AGC;  Oppenheim  has  applied 

generalized  linear  filtering  techniques  to  automatic 

correction  of  dynamic  range  limitations  in  photographic 
(43) 

images ;v  and  Fry,  et  al.  have  applied  computer 
control  to  the  problem  of  gain  adjustment. 

The  AGC  objectives  being  sought  for  in  sector¬ 
scanning  echocardiography  are  grey  scale  imaging  display 
and  avoidance  of  image  artifacts  caused  by  the  AGC. 
Optional  ability  to  reproduce  the  full  dynamic  range 
when  desired  has  been  implemented  because  of  research 
objectives.  Digital  signal  processing  is  used  to  permit 
the  study  of  different  processing  and  display  techniques. 
This  AGC  approach  consists  of  consecutive  functions: 
reversible  signal  compression  to  permit  A/D  signal 
conversion;  and  signal  compression  to  satisfy  the 
restricted  limits  of  display  dynamic  range. 

In  this  approach,  the  receiver  digitally  derives 
received  signal  power  statistics  to  control  the  TVG 
amplifier,  as  shown  in  Figure  A  .1.  The  purpose  of  the 
TVG  is  to  have  constant  A/D  input  level  by  using  known 
but  adaptive  gain.  The  TVG  is  a  piecewise-linear 
continuous  dB  curve  which  is  commanded  by  a  digital 
processor  in  terms  of  an  initial  value  and  15  slopes. 

The  maximum  range  is  typically  15  centimeters;  hence,  a 
different  slope  can,  in  this  case,  be  specified  for  each 
centimeter  increment  of  range.  The  data  needed  to  derive 


Figure  A.l  Adaptive-TVG  Pulse  Echo  Receiver 
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these  TVG  parameters  is  taken  during  a  set-up  scan  prior 
to  the  scans  for  diagnostic  images. 

A  window  is  used  to  cover  an  increment  of  range 
(and  also  of  angle  in  the  scanning  echocardiographic 
systems) .  Such  a  window  finds  use  both  to  calculate 
the  TVG  gain  profile  and  to  compress  the  dynamic  range  of 
the  image  for  display.  Figure  A. 2  illustrates  a  two- 
dimensional  image  raster  that  is  typical  of  a  scanning 
echocardiograph.  The  figure  includes  a  window  which  in 
size  is  9  degrees  by  m  millimeters.  To  calculate  the  TVG 
profile  for  a  two-dimensional  image,  one  can  make  a  choice 
between  different  window  shapes  and  signal  statistics. 

The  selected  statistic  in  dB  units  then  is  used  for  a 
least-mean-square  curve  fit  to  a  piecewise-linear  curve. 

The  piecewise-linear  curve  is  readily  converted  into  a 
TVG  command  function  which  controls  the  TVG  receiver  gain 
as  3hown  in  Figure  A.l.  The  steps  followed  in  TVG  profile 
generation  are; 

1.  Acquire  data  with  initial  gain  profile. 

2.  Combine  gain  and  receiver  output  to  yield 
echo  data. 

3.  Perform  weighted  LSF  of  logarithmic  echo 
data . 

4.  Convert  LSF  parameters  of  gain  profile. 

5.  Acquire  data  with  new  gain  profile. 

After  the  echo  signal  from  the  TVG  amplifier  output 


is  converted  into  digital  form,  the  signal  can  be  merged 
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with  the  stored  TVG  gain  commands.  This  restores  the 
signal  samples  to  the  initial  dynamic  range,  unaltered 
by  the  TVG  gain  profile.  In  an  echocardiogram,  almost 
all  look  directions  of  interest  contain  a  region  of 
blood,  so  that  some  small  duration  of  signal  contains 
only  backscattering  from  blood.  The  smallest  value  of 
the  selected  signal  statistic  (e.g.,  noncausal  r.m.s. 
value  taken  over  the  windowed  aperture)  will  occur  at 
such  a  region.  A  proportional  signal  threshold  can  then 
be  defined  at  this  range;  the  threshold  amplitude,  in 
general,  is  then  scaled  according  to  the  square  of  range, 
corresponding  to  the  divergence  of  sound  waves.  This 
adjusted  threshold  technique  has  been  applied  to  a  one¬ 
dimensional  (range  versus  amplitude)  segment  of  video 
heart-echoes  in  Figure  A. 3.  Here  the  dynamic  range  of 
the  signals  has  been  restored,  and  crosses  mark  the  peak 
echo-values  which  exceed  the  threshold.  Processing  for 
display  can  be  seen  from  Figure  A. 3  to  consist  of  several 
steps;  signal  detection  where  echoes  of  significant 
amplitude  are  identified,  clustering  where  separate 
tissue  structures  are  identified,  and  amplitude  compression 
where  each  separate  tissue  structure  is  scaled  so  that 
echoes  of  different  amplitudes  within  the  structure  are 
assigned  to  appropriate  brightness  levels,  but  all 
amplitudes  are  visually  discernible.  Figure  A. 3  gener¬ 
alizes  to  a  contour  map  in  the  case  of  two-dimensional 
scanning. 


The  application  of  adaptive  TVG  to  improve  the 
accuracy  of  the  A/D  converter  will  be  illustrated,  and 
an  example  will  be  given  of  an  amplitude  compression 
algorithm  for  image  display.  These  examples  will  use 
snapshots  of  A-mode  (amplitude  versus  range)  echocardio- 
graphic  data.  Generalizations  to  scanning  images  are 
straightforward.  These  demonstration  data  were  taken  by 
recording  the  transducer  output  of  a  commercial  M-mode 
echocardiographic  system.  The  test  set-up  is  shown  in 
Figure  A. 4.  The  test  data  exclude  approximately  the 
first  2.5  cm  of  range;  this  does  not  restrict  the 
validity  of  the  TVG  and  compression  algorithm  results. 

This  has  been  checked  in  subsequent  tests.  The  TVG 
amplifier  is  able  to  adjust  over  a  gain-range  greater 
than  90  dB. 

A  TVG  adaptation  algorithm  is  illustrated  in 
Figure  A. 5.  The  merged  TVG  gain  and  amplifier  output  is 
represented  in  the  upper  left-hand  curve.  This  represents 
the  envelope  of  the  RF  signal  at  the  TVG  amplifier  input. 
The  sample  signal  vector  in  dB  is  designated  by  y.  The 
diagonal  weighting  matrix  W  selects  the  peak  values  of  y. 
The  graph  of  y  and  the  best  mean-square  curve  fit  of  a 
piece-wise  linear  graph  to  Wy  are  shown  in  the  lower  left 
of  Figure  A. 5.  The  derived  gain  profile  is  shown  in  the 
upper  right,  and  the  lower  right  shows  the  TVG  amplifier 
output  which  results  from  using  this  gain  profile.  The 


Figure  A. 4  Echo-Level  Recording  System 


Figure  A. 5  Generation  of  TVG  Gain  Profiles 


output  has  been  rendered  relatively  stationary  compared 
with  the  input.  The  equation  representing  the  solution 
for  the  parameters  a  of  the  mean-square  fit  is 
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a  =  (FTWF)_1  FW  y  , 

where  F  is  the  matrix  such  that  the  sample  vector  of  a 
piece-wise  linear  curve  y  is: 

y  =  Fa 

for  a  given  vector  a 

Figure  A. 6  is  a  flow  diagram  of  an  amplitude 
compression  algorithm  used  for  echo  preparation  prior 
to  image  display.  The  application  of  the  algorithm  is 
shown  for  three  transducer  beam  directions  in  Figures  A. 7 
through  A. 9.  The  window  size,  (2N+1)  signal  samples,  is 
5.6  mm  in  these  examples.  The  window  in  this  algorithm 
serves  to  cluster  echoes  from  the  same  structure,  and  the 
r.m.s  values  taken  over  the  window  form  the  basis  for 
scaling  the  amplitudes  to  fit  within  the  display  bright¬ 
ness  range. 

Other  amplitude  compression  and  enhancement  methods 
can  be  used  with  this  system  implementation.  For  example, 
the  clustering  illustration  of  Figure  A .3  can  be 
implemented  by  thresholding  on  amplitude  and  identifying 


Figure  A. 6  Display  Algorithm 
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Figure  A. 9  Three-Dimensional  Case  of  Window-Normalizing  Algorithm 


peaks,  then  thresholding  successive  range  intervals  to 
identify  clusters,  followed  by  comparing  successive 
amplitudes  within  a  cluster  to  identify  regions  for 
amplitude,  or  brightness,  adjustment. 

The  design  of  a  system  for  clinical  use  will 
combine  the  two  adaptive  steps  of  TVG  signal  compression 
and  image  grey-scale  adjustment  into  one  real-time  TVG 
and  thresholding  operation  after  the  set-up  runs.  This 
will  be  done  within  the  system  organization  of  Figure  A.l 
by  design  changes  of  the  TVG  amplifier.  The  result  will 
be  processing  with  display  available  to  the  operator  with 
minimum  time  delay. 
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